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IV. On an inequality of long period in the motions of the Earth and Venus. 
By George Biddell Airy, A.M., F.R. Ast. Soc, F.G.S., late Fellow of 
Trinity College, Cambridge, and Plumian Professor of Astronomy and Ex- 
perimental Philosophy in the University of Cambridge. Communicated by 
Sir J. F. W. Herschel, F.R.S. 8fc. 85c. %c. 

Read November 24, 1831- 

J.N a paper " On the corrections of the elements of Deiambre's Solar Tables," 
published in the Philosophical Transactions for 1828, 1 stated that the compa- 
rison of the corrections in the epochs of the sun and the sun's perigee given by 
late observations, with the corrections given by the observations of the last 
century, appeared to indicate the existence of some inequality not included in 
the arguments of those Tables. As soon as I had convinced myself of the ne- 
cessity of seeking for some inequality of long period, I commenced an exami- 
nation of the mean motions of the planets, with the view of finding one whose 
ratio to the mean motion of the earth could be expressed very nearly by a pro- 
portion whose terms were small : and I did not long seek in vain. 

It is well known that the appearances of Venus recur in very nearly the same 
order every eight years : and therefore some multiple of the periodic time of 
Venus is nearly equal to eight years. It is easily seen that this multiple is 
thirteen : and consequently eight times the mean motion of Venus is nearly 
equal to thirteen times the mean motion of the Earth. According to Laplace, 
(Mec. Cel. liv. vi. chap. 6.) the mean annual motion of Venus is 650S-198; 
that of the Earth 3998*993. Hence 

8 X mean annual motion of Venus. . . . = 52018-584 
13 X mean annual motion of the Earth . . =5199 -909 



Difference = 1 *675 

The difference is about ^ of the mean annual motion of the Earth ; and it 

k2 
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implies the existence of an inequality whose period is about 240 years. No 
term has yet been calculated whose period is so long with respect to the 
periodic time of the planets disturbed *. The probability that there would be 
found some sensible irregularity depending on this term, may be estimated 
from this consideration ; that in integrating the differential equations, this 
term receives a multiplier of 3 X 13 X (240) 2 , or about 2,200,000. 

On the other hand, the coefficient of this term is of the fifth order (with re- 
gard to the excentricities and inclinations of the orbits). The excentricities of 
both orbits are small. And it is remarkable that in the present position of the 
perihelia, the terms which would otherwise produce a large inequality destroy 
each other almost exactly. The inclination however is not so small ; and upon 
this the existing inequality depends principally for its magnitude. 

The value of the principal term, calculated from the theory, I gave in a post- 
script to the paper above cited. I propose in the present memoir to give an 
account of the method of calculation, and to include other terms which are 
necessarily connected with the principal inequality. 

PART I. 

PERTURBATION OF THE EARTH'S LONGITUDE AND RADIUS VECTOR. 

Section 1. 
Method adopted for this investigation. 

1. The motion of a disturbed planet may be represented by supposing it to 
move, according to the laws of undisturbed motion, in an ellipse whose dimen- 
sions, and position are continually changing : the epoch of the planet's mean 
longitude at the origin of the time being also supposed to change. Putting a 
for the semi-axis major ; e for the excentricity ; zs for the longitude of perihe- 
lion ; n for the mean motion in longitude in a unit of time ; s for the epoch, or 
the mean longitude when t = ; (all which are variable) : m for the mass of 
the planet (Venus) ; f£ for the sum of the masses of the sun and planet ; and the 
same letters with accents for the same quantities relative to another planet (the 

* The period of the long inequality of Saturn is only about thirty times as great as the periodic time 
of Saturn. 
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Earth) ; the variation of the elements of the second planet's orbit will be given 
by the following equations : 

da' __ 2 ri d* dR 
dt ~ [>> de' 

dn' ___ 3 n B a> dR, 
dt — + pf "d7 

d^L _ ^JL n _ <2\ H*L _i_ n ' a '( l —**)* (dm , dR\ 
dt ~ ~ p!e> ^ 6 ' dJ "+" j»V W + d«7 

dvt\ __ n' a' . lQ> i dR 

dt ~ ~~ JJ7 \ l ~ e ) "de 1 

de' 3n'*a> dR . , In' a n dR n' a' ( ., ,„i „ ,„\ dR 

Ji = - ~^r •!!?•*+ ~y-'d^-V7\S 1 -^)- ( l - e2 ))d? 

m (a/ x + y'y) m . , , . 

where R or v+ir+*f ~ vw^wTw^wrv^m 1S ex P anded in 

terms depending on the mean motions of the two planets. These expressions 
are true only on the supposition that the actual orbit of m! is in the plane of 
x y, or is so little inclined that the square of the inclination may be neglected. 
The values of a', e', &c. on the right-hand side of the equations ought in strict- 
ness to be the true variable values. But it will in general be sufficiently accu- 
rate to put for e' the value E which it had near the time for which the investi- 
gation is made, and to consider it as constant: or at any rate the expression 
E + F t, where F is the mean value of its increase when t = : and similarly 

da' de' 
for the others. Determining thus the values of -^, -jj, &c. and from them 

those of d, e', &c, they are to be substituted in the expressions 
r' = a' 1 1 + \ e' 2 + (-e'+ -| e» - &c.) cos (n't + 4 - v') 

+ (~~ ¥ ^ + &c - ) cos (2 «'* + 2 «* — 2 O + &c. j 

v ' — n't + f!+(2e'-^ e' 3 + &c.) sin (n' t + s' - «r') 

+ (|- e' 2 - ~ e' 4 + &c.) sin (2 n't + 2 4 - 2 w 7 ) + &c. 
and the true values of the radius vector and longitude are obtained. 
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2. When (as in the present instance) the inequality is so small that we may 
be satisfied with the principal part of it, we may in the expressions omit the 
powers of e'. Thus we have 



da' 2n'a'* dR 
dt ~ fjj ' de' 










dn> 3n®a' dR 
dt ~" + ft' ' de> 


















de' _ ticl dR 

dt ~ ' /x'e' ' d^ 










dm' n'a! dR 
dt ~ ~ p'ef ' le 1 










de! Sn'^a! dR 2n!a'* 
dt —~~ fj ' de' • f + /*' 


dR 

' da' 


1 
2 


n' a' d 


dR 


V-' 


' de 1 



3. Hitherto this method has been actually used (I believe) only for the cal- 
culation of secular variations. But it can be applied with great advantage in 
almost every case : and in the instance before us it is particularly convenient, 
as it requires only the development of a single term. For if in the development 
of R we take the terms depending on cos { 13 n't — 8 n t + A}, whose coefficient 

d cv d ft d s 

is of the 5th order, it will be found that -rr, -jj, and jj, are of the 5th order, 

de' dts' 

jj of the 4th order, and -jj- of the 3rd order. Integrating these expressions, 

and substituting them in the formula for v', there will be produced terms of the 

forms ^joJzsnf sm {13ra'£— 8nt+B} and 13 n , q _ Q - sin { 12 n't — 8 n t + C}, 

where p is of the 5th, and q of the 4th order. And a little examination will 

show that no other argument will produce terms of the same or of a lower 

oi'der, which are divided by the small quantity 13 »' — 8n : inasmuch as this 

da' 
divisor is introduced only by integration of the expressions for -j7> & c - O ur 

object then at present is to select in the development of R all the terms of the 
form A cos {13 n' t — 8nt-{-B}. And as the inequality which we are seeking 
will probably be small, We may confine ourselves to those terms in which the 
order of the coefficient is the lowest possible : that is, to terms of the 5th order. 
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Section 2. 

On the abridgement which the development admits of, and the notation which it 

permits us to use. 

4. Let 6 be the longitude of the node of the orbit of m (Venus), and <p its 
inclination : the orbit of m' (the Earth) being supposed to coincide with the 
plane of xy. Let v, the longitude of m, be measured* by adding the angular 
distance of m from its node to the longitude of the node. Then v — d is the 
distance of m from the node. Let r be the true radius vector of m : then 

a?' = r'. cos v' 

y = r'. sin v' 

x = r {cos (v — 6) . cos 6 — sin (v — 6) . sin d . cos <p\ 

y =r {cos (v — 6) . sin 6 + sin {v — 6) . cos 6 . cos <p} 

% = r . sin (v — 6) . sin p 

Substituting these, the expression for R becomes 



mr' 



- j cos (vf — 6) . cos (v — 6) + cos <p . sin (w' — 4) . sin (v—d) > 



m 



a/ J /* — 2 r'r (cos (i/ — 0) . cos (v — &) + cos p . sin (»' — 0) . sin (v — S)J + r~ I 

in which it must be remarked that r and v, when expressed in terms of t, will 
not involve the constants 6 and <p. This may be changed into 

l~ < cos (v'~ v) — sin 2 -|- . cos («'— w) + sin 2 -| cos («' + w — 2 4) j- 



r 



?ra 



,/ J r J% — 2/r.cos(z/— w)+r 3 + 2r'r.sin 2 -^ cos(z/-u)— 2r , r.sin 2 -|-cos(i/ + a— 20) i 
or, 

— r COS (V — V_) — ^rpg _ Sr> ; — ^ - ^ - pj 

<b, f mr'r mr' 1 

+ sin 2 -|- {cos (v'— v)— cos (v'+v-™2 0} . -{ 7— ; ~ — ,.» -r > 

^ 2 ■ v ' v l{r' 8 -£r'r.cos(*/-7>)+r 3 p r » J 

* w, the longitude of the perihelion of m, must be measured in the same manner. 
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mr 1 r 



-Q- sin 4 —■ < cos (v' — v) — cos (v'r\r v — 26) > 



{r>* — % r 'r . cos (%/—v) + r*Y 

+ &C. 

5. The first line of this maybe expanded in the form 

f 1 (0) (\) , (2) , 1 

- m J ^ T\ + T^ cos (v'-v) +r^ cos (2 «' - 2 1>>+ &c. j- 

(o) (i) • 

where r. , I\ , &c, are functions of r' and r. We must then express r 1 and 

(0) (1) (2) 

r in terms of nt and rc *, and must substitute these values m I\ , I\ , I\ , &c. 

and must express v' and v in terms of n't and n tf ; and on multiplying the re- 
spective expressions we shall have the development necessary for our method. 

6. Now upon expressing" r in terms of n't, the following remarkable law 
always holds : The index of the term of lowest order in the coefficient of such 
an argument as cos (p n't + A), is p. The same is true with regard to the de- 
velopment of r, v', and v. 

7- Now such a term as A cos {I3n't— 8nt+B} can be produced only by 

the multiplication of C P^ fkn't — knt-\- k&'—ks\ , (from the first term in the 
development of cos (kv' — kv)\ with ^ (13 cv> ft) (n't + s' — -as') and 

sin (S^k) (nt -\- & — ■m) /'occurring in the development of k v' — k v, or of T T Y 

The largest term in the coefficient, according to the rule just explained, will 
be of the order whose index is the sum of 13 cv> k and 8 cv k. Now if k be < 8, 
as for instance if k be 7, the index of the order is 6 + 1 = 7, or the term is of 
the 7th order, and therefore is to be rejected. And if k be > 13, as for instance 
if k = 14, the index of the order is 1 + 6 = 7, and the term is to be rejected. 
But if k be 8, or 13, or any number between them, as for instance 10, then 
the order of the term is 3 + 2 = 5, and the term is to be kept. It appears 
therefore that the only terms which we shall have occasion to develope, are 

(8) (9) (13) 

V i . cos (8 v — 8 v), r t . cos (9 tr — 9 v), &c. as far as T } . cos (13 v' — 13f) 
inclusively. 
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8. Supposing then k to be not less than 8 nor greater than 13, the term 
(kn't — knt+ke'-ks\ must be multiplied by ™* ((13 - k) (n't + s' — z?') 

+ (k — 8) (nt + i — ■ut)\ in order to produce a term of the form A cos (13 n't 
— 8 n t -\- B) whose coefficient is of the 5 th order. The latter factor must have 



cos 

sili 



arisen from the product of two such terms as e' ' . • (13 — k) {n't + s' — ©•') 



and e .™* (k — 8) (n t + g — w). The expansion of such a product will always 

produce two terms, one of which has for argument the sum of the arguments 
of the factors, and the other has the difference of the same arguments. The 
point to which I wish particularly to call the attention of the reader is this : 
The term of the product depending on the sum of the arguments is the only 
one which is useful to us. For instance ; the product of e' 2 . sin 2 (n't + s' — d) 
ande 3 .sin3 (nt + s-vr) will be — J^^.cos (2n'l + 3nt-{-2e' + 3s — 2zs'— 3m) 
-f \ e' 2 <? 3 . cos (2 n't — 3 n t + 2 s' - 3 s — 2 rar' -j- 3 w) ; the combination of the first 
term with cos (\1 n't— lint -\- lis'— lis) will produce a term of the form 
A cos (13 n't — 8nt-{-B) whose coefficient is of the 5th order: the second 
term will not produce a term of that form. We might choose terms, as 
e' . sin (n't + s' — «r') and e 6 . sin 6 (n t + e — z*) such that the part of the product 
depending on the difference of the arguments, or ^ e e c > . cos (n't —6nt -\-s'—6s 
— ot'+6ot) combining with such a term as cos (14 n't— 14 n t -j- 14/— 14 s), 
would produce a term of the form required : but its coefficient would not be 
of the 5th order. It is equally necessary to remark that, in multiplying the 

term thus selected by C P^ (k n't — k n t + k s' — k s\, we again preserve only that 

part of the product depending on the sum of the arguments. 

9. On the circumstance that, in taking the product of two circular functions, 
we have to retain only the term whose argument is the sum of the arguments, 
depends the principle of our notation. For whenever (in an advanced stage 

of the operations) such a term as ^ ( 2 n't -\- 3 n t + 2 i + 3 s — 2 w' — 3 sA oc- 
curs, we shall know that, being formed in accordance with this rule, it must 

MDCCCXXXII. h 
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have arisen from the product of e' 2 J? s (W/+2£ , -2sr')and e 3 ^ (3«*+3s— 3cr ) ; 

its coefficient therefore can only be e' 2 e 3 . And conversely, from seeing this 
coefficient, we should be certain that the argument would be 2 (n't -\-4 — w') 
+ 3 (rc * + * — OT )- Instead therefore of writing 

e' 2 * 3 . cos (2 »'* + 3 n t + 2 s' + 3 s — 2 w'— 3 sr) 

we might simply write 

e' 2 e 3 . cos 

omitting the argument entirely. But it will be found more convenient to re- 
tain the figures in the argument, writing it thus, 

e' 2 e 3 .cos(2 + 3) 

the first figure being always appropriated to the accented argument. And 
when this term is multiplied by cos (1 1 n't — 1 1 n t + 1 1 4 — 1 1 4) or cos ( 1 1 — 1 1 ), 
we may write down the result 

1^^.008(13-8) 

without any fear of mistake. For we know that the argument must have been 
produced by adding 2 {n't -j- 4 — to-'), 3 (n t + s — ™), and 1 1 (n't —nt-\-4 — s), 
and thus when a result is obtained the term can be filled up. 

10. If we examine the second line in the last expression of (4), it is easily seen 

that sin 2 ~, a quantity of the second order (considering sin -x- as of the same 

order with d and e) enters as multiplier into two terms : of which the first, or 

sin 2 -|- . cos (v' — v), when developed will have in every term one part of the 

argument produced by a subtraction ; and therefore, when combined with the 
expansion of the term multiplying it, will produce terms cos (13 — 8) of the 
7th order at lowest ; the first term therefore is useless. But the second, or 

— sin 2 -|- . cos (y' + v — 2 6), is exactly analogous to e 2 cos (v' + v — 2 is), which 

would arise from the product of e 2 cos (2v — 2 w) and cos (v' — v), and to which 
all the preceding remarks would apply; and examination would show that in 

the development of this term, in which products of sin 2 -£ with powers of e' and e 
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will occur, the same rule must be folio wed, namely, that the only useful terms 
in the products are those in which the arguments are added. And whenever 

sin 2 -|- occurs in the coefficient, — 2 6 occurs in the argument ; so that there 

will be no possibility of mistake in using the notation described in (9). 

11. On examining the third line in the last expression of (4), it will be 

seen in the same manner that the only part of ■— — sin* -|- < cos (v' — v) 
— cos (if + v — 2 6) > to be preserved is — — sin 4 -| . cos (2 v' -f- 2 v — 4 6). 

The same remarks apply to this term as to the last ; and for a similar reason 
the notation of (9) may be used without fear of mistake. 

12. By the use of this notation we may in some instances materially shorten 
our expressions. For instance, we might have the terms 

Fe'e 4 .cos (n't + Ant + s' + 4s — «/ — 4w) 
+G<?Vsin 2 -1 . cos (n't + Ant + { + 4e- w'-2 w — 2 6) 

+ H d sin 4 1- . cos (n't + 4 n t + { + A s — z*' — A 6) 

All this would be expressed without the possibility of mistake by the follow- 
ing term, 

(F 4 <* + G J ei sin 2 ± + H e' sin 4 -IV cos (1 + 4) . 

The utility of such abridgments, and the quantity of disgusting labour which 
they spare, can be conceived only by those who have gone through the drudgery 
of performing the actual operation. 

13. It is only necessary to add that when we have, for the coefficient of a 

cosine or sine, a series proceeding by powers of e', e, sin 2 -^ , &c. we may always 

neglect all after the lowest power. For instance, the correct expression for 

v is 

nt-\-t 

(1 5 \ 

2e— -r^+oj^ 5- & c -) sin («£ + «— zs) 

l2 
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+ (LeL— —e i +8>cc.\sm(2nt + 2z — 2sr) 

+ (l|«3-||e 5 +&c.)siii(3n/ + 3s-3TO-) 
+ OS- e 4 -&c.) sin (4 n * + 4 s - 4 w) 

+ (-^ e 5 -&cA sin (5 ra* + 5 s-5 or) 

+ &c 
but for our purposes it will be sufficient to take v = (0 + 1) + 2 e . sin (0 + 1) 

+ !-«*. sin (0+2)+ ^.8111(0 + 3)+ ^.sin(0+4) + |g>.sm(0+5). 

For none of the terms can be of any use to us till they are multiplied, so that 
the largest term of the coefficient is of the 5th order ; and then all the other 
parts will be of a higher order. 

14. Putting/ for sin -§-, it will be seen that (in conformity with the remarks 

4* 

in this section), the terms of R to be developed are 



"/ 



\/{r' % — Qr'r. cos (w' — v) + r 2 } 
mr'r . cos (rf + v — 2 0) 



{r' 2 — 2 / /• . cos (t/ — v) + r 3 }" 2 



3 fi ffl / 2 r a . cos (2 1>' + 2 w — 4 fl) 

"""" 2jT/" • ~ — — — _ 

{r H — 2 r" r . cos («/ — v) + r 9 } T 

Section 3. 
Expansion of cos {kv' — kv), to thejifth order. 
15. By (13) the value of kv' — kv is 
(A -A) 
+ 2&e'.sin(l+0)-2£e.sin(0+l) . (A) 

+ ~ke , Ksm(2 + 0)-^-ke' i .sm{0 + 2) • ( B ) 
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+ ^ke'K sin (3+0) - 1|^. sin(0+3) , .(C) 

+ ^?Ae , *.sin(4 + 0)-™*A**.sm(0+4) (D) 

+ ^A^sin(5 + 0)-^gA^.sin(0 + 5) (E) 

The cosine is 

cos(A-A).cos(A+B + C + D + E) 

-sin(A-A).sin(A+B + C + D + E) 
or 

cos (A- A). 1 1 g + ^ j 

• ^7 i\ lAJ.RjLriT^_T7 A 3 +3A 2 B + 3A a C + 3AB a , A 5 ") 
— sm(A— A) „ jA+B + C + D + E g + I20f 

omitting all products of an order above the fifth. 

16. In expanding the powers of A, B, &c, and in multiplying the expan- 
sions by cos (k—k) and sin (A — A), the rules of (8) must be strictly followed. 
Thus we find at length for the value of cos (kv' — k v) : 

Principal term, 
cos (A — A) 

Terms of the first order, 
+ Ae'. cos(>fc+l— A)— A e. cos (A— A— 1) 

Terms of the second order, 

(y k2 + i k )^ K cos (^+2 - A) - A 2 e' e . cos (a+1 - A^l) 

+ (t k% ~ i k ) e2 - cos (A- A^2) 

Terms of the third order, 
(lA3+4A 2 +^A) e '3. C os(Aq^-A) + (-^A3_4A2) e ^.cos(A+2-A=~ 1 ) 
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+ i~ &- -| A 2 ) e 1 e* . cos (a+1 - k^2\ 

+ (..^*l + i.A._.g*) a 8. C 0B(*-J3s) 

Terms of the fourth order, 
+ (- ^W-^tf-£ 4 k*)e*e.co S (k + 3-J^l) 

+ (■? **- M ^V 2 e2 • cos (^+^ ~ ^~^) 

Terms ofthejifth order, 
+ (f 2 A5 +4 A4 -S^-S^)^^-^s(FF3-FT2) 

+ (s k5 - f 6 **+ i->- ii *) ^ • cos (*+* - ^) 

This development includes every argument whose coefficient is of an order 
not exceeding the fifth, The coefficients however here exhibited are only the 
first terms of the series which represent the complete coefficients. 
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Section 4. 

(k) 
Expansion of •— I\ , to the fifth order. 

17- We suppose — ^^ _ g r , - - cos (p , _ ^ - r ^ , the first term in the expres- 
sion of (14), to be expanded in the form 

— -^m,r i —mT i . cos (v'—v) — mT i . cos (2 v' — 2 v) — &c. 
-mr, . cos(£t>' — &v) — &c. 
where I\ , 1^ , &c. are functions of r 1 and r only. 



Let 



*/{a m — la' a cos (t/ — v)+ a 8 } 

= 4 cf + cf .cos(v'-v)+cf .cos{2v'-2v)+&c. + cf .cos(kv'-kv)+kc. 

then I\ is the same function of r 1 and r that C^ is of a! and a. Consequently, 
if r' = a' (1 + q'), r = a(l-\-q): and if for convenience we use the notation 

(m, n) CL 
to express that which is commonly written 

-,'m 



d m . a n 



da m .da n 
{k) 
we shall have for — T i the following expression : 

- <<>,< . ,- (,,,)cf . rt- (2,< . £- (3, 1)C f . ^ - (4,1)4". & 

- (0 ;2 ) C f . t _ (1;3)c f . if. (2;2) cf . «_ (3;2)c f . « 

-<0,3)cf.£ -( 1 ,3)cf.^-(2,3)cf&' 

-tM(f.fi -(.,4)cf.^ 
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18. The value of r, contracted according to the system of (13), is 
ai 1 -e. cos (0+ 1) — \ e 2 . cos (0 + 2) — -f-e^cos (0 + 3) — ^-e 4 . cos (0 + 4) 



e 5 . cos 



(0 + 5)} 
whence q = 



125 
384 



- e cos (0 + 1)— y e 2 . cos (0 + 2) — -g- e 3 . cos (0 + 3) — \ e* . cos (0 + 4) 
-384^. COS (0 + 5) 

and a similar expression holds for q'. Substituting these in the expression 
above, and following strictly the precept of (8), we find for the development 

of-4 ; , 

Principal term, 

- c'f 

Terms of the first order, 

(k) (k) 

+ (1,0) C, . J ccs (1 + 0) + (0,1) Cj . e cos (0 + 1) 
Terms of the second order*, 
(|(1,0)- ~ (2,0) j cf .^ cos (2 + 0)-- i-(l,l) cf . e'ecos (1 + 1) 

+ ( ± (0,1) - i- (0,2) } Cf . e 2 cos (0 + 2) 

Terms of the third order, 
(4(1,0) -4 (2,0) + 55(8,0)}cf.^coB(3+O) 
+ (- |(1,1) + -|-(2,l)}cf.^ e cos(2+l) 

* In tHs and the succeeding expressions, when a cosine is multiplied by the sum of several diffe- 
rential coefficients of C^ , the symbols of differentiation are bracketed together, and C^ is put at the 
end of the bracket. 
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+ {- i-(l,l) +| (1,2) jcf.e'e* cos (1 + 2) 
+ (4(0,1) -4(0,2) + ^(0,3)jcf.^cos(0 + 3) 

Terms of the fourth order, 
{¥( 1 ' )-l( 2 ' ) + ^( 3 ' )-i^( 4 '°)} c f- e ' 4cos ( 4 + ) 

+ {--S( i ' i )+4( 2 ' i )-i( 3 3 i )} c f- e ' 3ecos ( 3 + 1 ) 

+ {-4(l,l)+^(2,l)+^(l,2)-^(2,2)}cf. e '^2 C0S (2 + 2) 
+ {-^(l,l)+4(l J 2)-^(l,3)}cf. e Vcos(l+3) 

+ (1 (°>V ~ j(°> 2 ) + ^(°' 3 ) " l^(°^)} C f • e * cos (° + 4 ) 

Terms of the fifth order, 

{ m 0»°) - I ( 2 >°) + & ( 3 >°) - gB ( 4 >°) + i^o ( 5 >°) } c f • *' 5 cos ( 5 + °) 

+ { - | (1,1) + 4 & 1 ) - 5 ( 3 > X ) + 381 (4.1) } C f • *' 4 ' C ° S ( 4 + 1) 

+ { - I (LI) + 16 ( 2 '1) + (I ( J ' 2 ) ~ W* OM) ~ 5§ ( 2 ' 2 ) 
+ ^(3,2)}cf.^cos(3 + 2) 

+ { ~ I (!'!) + £ ( 2 >*) + 16 0> 2 ) - 31 (^ 2 ) - 96 (M) 
+ i ^(2,3)}cf.«'Mcos(2 + 3) 

+ { - 4 (M) + T (^ 2 ) - k 0> 3 ) + 3§4 (M) } C f ■ ' •* c <* (1 + 4) 
+ { 384 P. 1 ) ~ 96 (°> 2 ) + 64 (°> 3 ) ~ 58"(0,4) + ^ (°> 5 ) j C i . * COS (0 + 5) 
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Every argument is included whose coefficient is of an order not superior to 
the fifth : but only the lowest order of each coefficient is taken. 



Section 5. 

Selection of the coefficients of cos (13 — 8) in the development of 

m 
~ a/ {V 2 — 2 r 1 r . cos (»' — v) + r 2 } 

19. For this purpose, as the general term in the expansion of 

77-5 — ■x—, t~i \~> — ¥, is —mTz . cos (k v — k v), we ought to mul- 

tiply together the expressions of (16) and (18), to multiply the product by m, 
and then giving different values to k to select those terms which have for argu- 
ment (13 — 8). But without going through this labour we may, when a value 
is assumed for k, select by the eye the terms required. As we have explained 
in (7), the values which it is proper to give to k are 8, 9, 10, 11, 12, 13. 

20. Thus we obtain the following coefficients of cos (13 — 8) : 

k=8. 



m 



f 639753 . . 178109 4217 , . 407 , . 

X ( 240" (°»°) + ~TW V 1 ' ) ~ "198 ( 2 '°) + 384 C 3 »°) 

-T9 L 2 ( 4 ' o )+34M)} c ^^• • .' (l (s) -*' 5 ) 



k = 9. 



m 



f 1955097 , , 88029 . , 217233 4041 

X ( 384 (°»°) ~ ~W ( lj0) + ~768~ CO* 1 ) + ^4^ ( 2 >°) 

- w 0' 1 ) ~ w < 3 ' ) + t§ (2 ^ + lis ( 4 '°) ~ <s (3*1) 

+ 7^( 4 ' 1 )} C i 9, ^' 4 « (L (9) .^e) 

(8) (8) 

* By (0,0) C is meant the same as C. . 
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k = 10. 

f 492625 , n „, , 86275 ., _ 53485 , n . 1925 . „. 
m >< { - -48" (°'°) + -W (^°) - "IS" (°^) - -IT ( 2 '°) 

, 9367 , , 2815 175 , 209 , . 493 
+ ~6T C 1 ' 1 ) ~~ ~W (°' 2 ) + "192 ( 3 >°) ~ ~32 V 2 ' 1 ) + "128 ( X > 2 ) 

+ iS & 1 ) " Bi ( 2 > 2 > + 381 ( 3 > 2 ) } 4 I0) • «'* # • • ' ( L(10) - e ' 3 e2 ) 

A=ll. 

_ f 498107 /rt . nN 20999 ,, n x ■ 52283 913 

m X J -4g— (0,0) - -g^- (1,0) + -g^- (0,1) + ^g (2,0) 

2231 « 97, . 115 M9 fn „, 

+ l2( 2 ' 2 )-S( 1 ' 3 ) + 3k^ 3 >} C * n) -« ,2 « 3 • • • (L (U) .e'2^) 

A = 12. 

f 20337 /A/ ,, , 6779 /1AX 211 7/„-i\ , 2119,, lX 
m X J J- (0,0) + -32- (1,0) ~ — (°,!) + -48~ t 1 ' 1 ) 

~ ^ (0,2) + ~ (1,2) " t (°, 3 ) + 61 (1' 3 > ~ S(M 

+ 4(l, 4 )}cf ; eV. . (L^.eV) 

A = 13. 

f240643. A .. , 24571 , AlN . 1219 „ _. , 235 
m X { "alo" (°' 0) + ~W~ (0' 1 ) + "4T (°'2) + 192 (0, 3 ) 

+ 3^4 («.4) + 3^5 (0,5)}cf. e s (L^.ea) 



The arguments of the cosines multiplied respectively by these coefficients, it 
must be recollected, are not similar. Their form will be determined by the 
considerations mentioned in (10). 

m2 
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21. The next term of R to be developed, by (14), is 



— m . -r ./ 2 . cos (v' + v — 2 $) 

{r 13 — 1 r'r . cos (»' — ») + t & Y' 

We shall put T 3 for the general term in the expansion 

fJ r (0) (1) (2) 

s „ - : : 7rL=i T i + r * .cos(t/-») + r/.cos(2t/-2»)+&c.; 

{r'* — 2r'r.cos(i/ — ») + r 4 }^ T T ^ 

And C 3 for the general term in the expansion 

— ■ a ^~ -_ = i C f + C ( 3 1} cos (v'-v) + C ( , 2) cos(2t/-.2t,) + &c. 

{a'2-Qa>a.cos(<if-v)+a?}- !r * » 

Section 6. 

Development off 1 . cos (if + » — 2 #), £o the fifth order. 

22. As the multiplier/ 2 is of the second order, we want cos (v' -\-v — 2 6) 
only to the third order. Now, by (13), if + v — 26 — 

(1 + 1) -26 

+ 2e'sin(l + 0) + 2<?.sin(0 + 1) (A) 

+ Xe' 2 .sin(2 + 0)+|e2.sin(0 + 2) (B) 

13 1 ^ 

+ - e '3sin(3+0) + I 2e3.sin(0 + 3) (C) 

Its cosine, as in (15), is 
(1 + 1-20). {l-^±|ABj_. sin(1 + 1 _ 2 , } | A + B + c _A3| 



cos 



Following the rule of (8) in the expansion, we find for the value of cos 
(v' + v — 2 6). 

Principal Term, 
cos (1 + 1 — 2 6) 
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Terms of the first order, 
+ e' . cos (2 + 1 — 2 6) + e . cos (1 + 2 — 2 6) 

Terms of the second order, 
+ -| e' 2 . cos (3 + 1 -2 4) + e'e.cos(2 + 2-20) + |-e 2 .cos(l +3-20) 

Terms of the third order, 
+ — e' 3 . cos (4 + 1 - 2 0) + |- e' 2 e . cos (3 + 2 - 2 0) 

+ |- e'e 2 . cos (2 + 3 - 2 6) + -^ e 3 . cos (1 + 4 - 2 6) 

On multiplying this by/ 2 it will readily be seen that/ 2 in the coefficient is 
always accompanied by — 2 6 in the argument, and that there is a necessary 
connexion between them. We may therefore omit 2 6 -, and thus we have for 
the development of/ 2 . cos (v' + v — 2 6) 

Term of the second order, 
/ 2 . cos (1 + 1). 

Terms of the third order, 
+ e'/ 2 . cos (2 + 1) + ep . cos (1 + 2). 

Terms of the fourth order, 
+ |- e' 2 / 2 . cos (3 + 1) + J ep . cos (2 + 2) + f e 2 / 2 . cos (1 + 3) 

Terms of the fifth order, 
+ T eV 2 • cos (4 + 1) + -| e' 2 e/ 2 . cos (3 + 2) + ~ e' e 2 / 2 . cos (2 + 3) 

+ i e 3/ 2 .cos(l+4) 
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Section 7. 
Development of cos (kv' — k v) .f 1 . cos (v' + v — 2 &), to the fifth order. 

23. We must multiply the expression in (16), (of which only the terms to 
the third order will be wanted), by the expression just formed, according to 
the rule of (8). Thus we obtain the following expression : 

Term of the second order, 



g-/*.cos(ft+ l-k-1). 

Terms of the third order, 
(|^+i)e7 2 .cos(yH^-^^)+(-Y^+|) e / 2c os(ATl-^ ; ^2). 

Terms of the fourth order, 
(t* 2 + TB^ + ^W 2 -cos(FR-£3T) 

+ (- \kz + 7 *-) e'e/ 2 . cos (A +2 — k - 2) 

+ (t* 8 - T6 k + To) e 2 / 2 .cos(I+l - £^3) 

Terms of the fifth order, 
(i^ + ^^ + l^ + §) e V 2 -cos(X+4-^l) 

+ (- i^ 3 - to ^ 2 + i * + to) e ' 2g / 2 - cos (^+3 - F=r 2) 

+ (t F - TO^ 2 ~ T* + to) e ' e2 / 2 • cos (^+2 - ^ r 3) 
+ (-^^ + TO fc2 -| fe + l)^/ 2 - C0S ( I + T -^ :ri ) 
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Section 8. 
Selection of the coefficients of cos (13 — 8) in the development of 

aJ At 

— m . -f 2 . cos (v' + v — 2 6). 

{r 12 — 2 r 3 r. cos (if — e) + > 8 } T 

(*) 
24. The general term of the expansion is — m . T 3 . cos (kv' — kv) .f 2 . cos 

"2" 

(t/ + w — 2 4). The expression for cos (k v — k v) .f 2 . cos (v' + w — 2 6) we 

CO 
have just found ; and the expression for — T, will be in all respects similar 

(k) (*) (4) 

to that for — I\ in (18), putting C 3 for C* . Observing that k cannot be 

less than 9 or greater than 12, and selecting for the different values of k the 
terms whose combination produces (13 — 8), we get the following coefficients : 

fc = 9. 



m 



X { ~ ^ (°>°) + ll (M>) ~ IB (2»0) + §S ( 3 '°) } C i 9) ' e V 2 • • (M ( V/*) 



fc= 10. 

» x {^(0,0) _ 52Z(i,o) + ^ ( o,i) + ^(2,0) - ^ (i,i) 

+ &<?,!)} C^.Pef* (M (10) .e'2e/2) 

& = 11. 

« X { - ~ (0,0) + ^ (1,0) - f (0,1) + i| (1,1) - 4 (0,2) 

+ 4(l,2)}c ( " ) . e 'e2/2 (M (U) .eV/2) 



& = 12. 



m 



X {i^ (0,0) + f (0,1) + 4 (0,2) + | (0,3)} Cf } . e3/2 . . (M^.eV 2 ) 
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The arguments of the cosines multiplied by these coefficients are not similar ; 
their forms may be found by the reasoning in (10). 
25. The next term of R to be developed, by (14), is 

— m.-r — r/ 4 • cos (2 v' 4- 2 v — 4 6). 

4 ^ 3 _ 2/ r .cos (»' -v)+ r«} x 

(*) 
We shall put r 5 for the general term in the expansion 

r ,2 r a (o) (i) , (2) 

. — — = k + r; . cos (»'-») +r; . cos(2t/-2»)+&c. 

{r't-QSr. cos (rf-v) + r*}* * * * V ' ' 

and C 5 for the general term in the expansion 

r<n„% (0) (1) (2) 

; = iC. + C. . cos (t/-w) + C. . cos (2 v'-2v) + &c. 



{a' 2 -2«'«.cos(»'-v) + « 2 }" ; 



Section 9. 

Development of cos (kv' — /c v) ./ 4 . cos (2 «' -f- 2 f — 4 0), to ^e j^ftA order. 

26. As the multiplier / 4 is of the fourth order, we need to develope 
cos (2 v' + 2 v — 4 4) only to the first order. Now by (13), 2v } + 2v — A6 = 

(2 + 2) — 4 4 

+ 4e'. sin (1 +0) + 4 e . sin (0 + 1) 

and consequently cos (2 v' + 2 w — 4 4) = 

cos (2 4- 2 - 4 4) — sin (2 4- 2 - 4 4) . {4 e' . sin ( 1 + 0) 4- 4 e . sin (0 + 1 ) } 

= cos (2 + 2 — 4 4) 

+ 2 e' cos (3 + 2 — 4 4) + 2 e . cos (2 + 3 - 4 4) 

Multiplying this by/ 4 it will be seen, as in (22), that we may omit 4 4 in the 
argument. Thus we have for the development of/ 4 . cos (2 v' + 2 w - 44), 

Term q/' the fourth order, 

/ 4 .cos(2 + 2). 
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Terms of the fifth order, 
+ 2 e'/* . cos (3 + 2) + 2 e/* . cos (2 + 3). 

27. This is now to be multiplied by cos (kv 1 — kv), the expansion of which 
has been performed in (16). Effecting this operation, we have for the deve- 
lopment of cos (k v' — kv) .f* . cos (2 v' + 2 v — A 6), 

Term of the fourth order, 
1 



2 



•/« . cos (A + 2 — A - 2) 



Terms of the fifth order, 
(y A + 1 ) e'/* . cos (FHK3 - F=2) +(--ij + l)e/4. cos (1+2 -1^3) 



Section 10. 
Selection of the coefficients of cos (13 — 8) m tfAe development of 

{^a _ 2 rf r . cos (t/ — »)■ + r 2 }* 



- w • T • ~ „ , ~ r — r7I ./*.cos(2i/+2w-4^). 



3 
28. We must suppose the expression of (27) to be multiplied by -jm, and by 

the expression for — T„ I which will be formed from that of (18), putting C 

for Cj J. Then giving to k different values, we must select the terms in the 
product whose argument is (13 — 8). It is easily seen that 10 and 11 are the 
only admissible values of k. Thus we get these coefficients; 

k—\0. 



m 



x{-|(0,0)+^(l,0)}cf ) .e'/4 (N £W) .e'/*) 



ft=ll. 



mx{f (0,0) + -~ (0,1) }c^\ef (N (U) .e/4) 
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29. The terms collected in (20), (24), and (28), form the complete coefficient 
of cos (13 — 8) in the development of R to the fifth order. The arguments of 
the cosines multiplied by the different series are all different ; so that there 

(8) 

are twelve different terms to be calculated. Using the symbols L , &c, the 
complete term is expressed thus : 

(8) 

L .e' 5 .cos {13 (n't + «') — 8(n* + g) - 5sr'} 

(9) 

+ L . e' 4 e . cos {13 (ri t + g') -8(»/ + t)-4»- «r} 

(10) 

+ L . e' 3 e 2 . cos {13 (n't + a') — 8 (n t + e) - 3 w' - 2 sr} 
+ L (H .^e 3 . cos {13 (n't + s') — 8 (nt + g) — 2w' - 3 w} 
+ L 12 . e' e 4 . cos {13 (n't + g') — 8 (nt + s) — w' — 4 w} 

(13) 

+ L . e 5 . cos { 13 (n' t + s') — 8 (w t + g) — 5 sr} 

(9) 

+ M . e' 3 / 2 -cos{13(rc'* + g') - 8(»*-f-g) -3w'-24} 

(10) 

+ M . e' 2 e/ 2 . cos {13 (ri t + e') - 8 (n t + 1) - 2™' - n - 2 0} 
+ M (U) . e' e 2 / 2 . cos {13 (n't + s') - 8 (n* + s) - «r' - 2t* - 2 0} 

(12) 

+ M . e 3 / 2 . cos {13 (rit + g') - 8 (nt + g) - 3 ™ - 2 0} 

(10) 

+ N . e'/ 4 . cos {13 (ri t + s') — 8 (nt + e) — v> - 4 6} 
+ N (U) . e/ 4 . cos {13 (rc'* + g')- 8(«*+g) -sr-40} 



Section 11. 

Considerations on the numerical calculation of the inequalities in the Earth's 

motion depending on this term. 

30. If we examine the expressions of (2), it will appear that the values of all 
may be deduced with little trouble from the terms above, except that depending 

on -jfl. Since a' enters only into the coefficients, j^ will be produced by 
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differentiating the coefficients and retaining the same cosines. The coefficients 

(8) 1 (8) (8) 

will be differentiated by changing (0,0) C 4 into — (1,0) C 4 , (3,2) Ci into 

1 (8) 3 (8) 

-j (4,2) C j + -j (3,2) C 4 , &c. Thus new terms will be introduced whose cal- 
culation is rather troublesome. It is desirable, then, to inquire whether it is 
probable that the term depending on j^, will be comparable in magnitude to 

the other term which has the same argument. 

31. Now if we put A.cos{13(rc'i? + a') -8(nt + i) + B] or Acos (13 — 8), 
for one of the terms, we find 

-^ = — 3 . 13 . — r . A . sin (13 — 8) 



whence 



»' = N' + (1 |^f|^,A. cos (13-8) 



(where N' is constant and = mean value of n') 

w = +3,13. -^7-A.#.sm(13-8)+-^-. jj cos (13 -8) 
whence 

, „, 3.13ra'*a' . , ,,„ „. , 3,13ra' 2 a' . . .,„ . 

g==E - (13^- 8re )^ A -^ COS ( 13 - 8 )+ (13 re '-8^ V A - Sin ( 13 - 8 ) 

+ ll3n'-8n)[J ' dd Sm V id ~ y J 
(where E' is constant and = mean value of s') 
and n't-\-s' (which, by (1), is the first term off') becomes 

N^ + E' + { (13 ^_ 8w)V , A+ (13?t> _ 8w) ^ .^jsin(13-8). 

The ratio of the two coefficients of the inequality sin (13 — 8) is 

39 n' f dA 

2 * I3n'-8n : a da! 



dA 
or nearly 4800 X A : a' j^. 



n2 
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It will be seen hereafter, that for any one of the terms whose union com- 
poses L , &c, d j^r is greater than — A, and that it may, on the mean of 

values, be said to differ little from — 1 2 A. This reduces the ratio of the terms 
to 400 : 1. Now though we cannot assert that the sum of one set of terms 
will have to the sum of the other set of terms a ratio at all similar to this, yet 
the great disproportion of the terms related to each other seems sufficiently to 

justify us in the a priori assertion that the terms depending on -r-> are not 

d R, 

worth calculating. It will readily be seen that the terms depending on -j-j are 

J T> 

still more insignificant than those depending on — r 

32. We stated in (1) that the variations of the elements would be sufficiently 
taken into account in the expression for R if we put E + F t for e, &c. ; which 
amounts to taking only the secular variations. There will be no difficulty in 
doing this for e 1 , e, «/, vr,f, and 6 : but if such terms existed in the approximate 

expressions for a 7 and a, they would require the use of the differentials -j^>~t~- 

But a 1 and a have no secular variations : and therefore these differentials are 
not wanted. We may therefore proceed at once with the numerical calcula- 

(8) (9) 

tion of the terms L , L , &c. 

Section 12. 

(0) (1) (2) (£) (i) (*) 

Numerical calculation of C , C , C , ofc, C, , C 5 , $c. to C . 

§ i i T T v 

33. If we put t — 2 w for v' — v, we have 

1 (°) (0 ( 2 ) 

- ; -75 , „ , — x ; — jt = i Ci — C x . COS 2 a 4- (X . COS 4 a> — &C. 

4/ {a'* + 2 a' a. cos 2 m + a 1 } 3 i i ' t 

Integrating both sides with respect to a, from a = to a = -^ and putting S^ 
for the symbol of integration with respect to a between these limits, 

<; I _JLr {0) 

°» * y'Ja' 8 +1a! a. cos «» + a 2 } ~~ 4 ^i 
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whence 

C 4 ~~ m S a • y{a' a + 2a'a.cos2co + a 3 } 
or, putting a for -^-, 

c. = — i • S 



J i ira' ' & « ,/{! + 2acos2a> + a 3 }" 



Now let sin */ = V{ t + aTco^ , + gj ; and «f = ~^f~ : after substitu- 

tion it is found that 

(0) 4 1 

C § = Za't 1 "*" "') * ^ * a/{1 + 2«'cos2«)'+a' 2 } 

_ . , . . „ sin 2 a/ „ 1 _ VT^aP' 

In the same manner, making sin *> = ^ {1 + 2rfcosft«/ + «'*} : a = i+^rT ?* : 

(0) 

and so on, we get for C 4 the expression 

The values of a', a", &c. decrease very rapidly; and when a (n) is insensible, 
S.w — tsA m fsr. becomes S„w . 1 or~. Consequently 

■/{ 1 + 2 « ( ; cos 2 ce w + aS ' } * 

4 0) = -|,(l+O (l+«") (l+«'") .&C. 

the factors being continued till a (n) becomes insensible. The calculation is very 

& &' (°) 

easy; for, if we make sin|3 = a, sinj8'= tan 2 ~, sin 0" = tan 2 -^, &c. then 1^ 

_ i. sec 2 L . sec2 £ . sec 2 ^ . &c. For Venus and the Earth (Mec. C61. liv. VI.) 

(°) l 
« or -7 = 0,7233323 : using this number in the calculation, C 4 = —, X 2,386375. 

34 Affain, , r s — ^— ; ^ — ; — a? = tCi. cos 2 <w — -x- C 4 (1 + cos 4 &») 

°*' n 6 ttlu ' \/{a' a + 2 a' a. cos 2 eo + a*} 2 4 2 4 v ' ' 

1 (2) 

+ — Ci (cos 2 a + cos 6 a) - &c. ; integrating between the same limits as 
before, 
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(i) _4_ „ cos 2 to 

W — — Vd ■ " ' V{ 1 + « « cos 2 eo + « 3 } 

Making the same substitution as in (33), there will be produced three terms ; 
of which one vanishes in the definite integral, the second is similar to the ex- 
pression of this article, and the third similar to that of (33). Making a simi- 
lar substitution in the second term, new terms are produced. Pursuing this 
method, it will be found that the only terms whose values are ultimately sen- 
sible are those which are similar to the expression of (33) : and at last we get 



-,(!) 



=4 0, -{ 



sin |8 sin |3 sin /3' sin /3 sin 0' sin /3" 



+ &c.j=-^,X 0,9424 137 



35. Putting % for v' — v, and differentiating with respect to % the logarithms 
of both sides of the equation 

1 1 (°) (0 (2) 

,; i s — rr- -, ' , 8 , =ttC, + Ci . cos r + Cx cos 2 y + &c. 

multiplying out the denominators, and comparing the coefficients of cos k%, 

(*+ 1) __ 2 k (\_ \ _(*) 2k- 1 (*- 1) 

H — 2k + 1 \ cc + a ) W ~ 21c + 1 C * 

where — -}-« = 2,1058226. Making A; successively 1, 2, 3, 4, &c, we get the 
following values : 



(o) 

C\ = 



c ( : } = 



X 2,3863750* 
X 0,9424137 



( 2 ) 

C, = 7 X 0,5275791 

(3) 

C^ = -j X 0,3233422 

( 4 ) 
C\ = -j X 0,2067875 

C^ =77 X 0,1355852 



(6) 1 

■C 4 =-,X 0,0903724 

(7) 1 

Cx = -f X 0,0609432 



cf = Jr X 0,0414571 

(9) 1 

C, = -r X 0,0283925 



(io) i 
C 4 = ? X 0,0195495 

Cj U) = -^ X 0,0135189 



(12) 1 

Cx =-7X0,0093812 



Cx 13 = ~ X 0,0065274 

(14) 1 

Cx = -r X 0,0045503 



(15) 1 

C i =^X 0,0031744 

(16) 1 

C 4 =-jX 0,0022123 

(17) 1 

Cj =-rX 0,0015356 



c 



(18) 



-r X 0,0010554 
a! 



* Laplace's numbers, which are somewhat diiFerent from these, are computed by the less accurate 
method of summing a slowly converging series. 
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(*) (*) 
36. For the calculation of the terms C 3 , C 5 , &c., we shall adopt the general 

notation 

_ ___ = - C, + C s cos x + C s cos 2 x + &c. 



, fa' „ . a ~) ' 

N /a'«.(--2cos % +^rj- 



which, it will be seen, includes those of (17), (21), and (25) ; and proceeding 
as in (35) we shall find this general equation 

And since 






N /« r «.{T- 2cos ^ + -FJ" 

(± + ._2cOB X )x— — p— J ,a1. + i 

X ' J°>a.{ T -2cosx + -a\ 

we find on substituting the expansions and comparing the coefficients of cos &%, 
C®-( 1 +«\C W _C < *' ,) -C ( * +1) 



(* + l) 
Removing C 4 + j by means of the relation just found (putting s -f- 1 for s) 

In nearly the same manner, 

(fc-1) 

Eliminating C s + j , 

r (t) - _ a (* + *-• i) i pft-Q ft-i «T + * r (*> 

■ +, ~ * ■(£-.)• * " * •(i-y * 

If in this we substitute the value of C s in terms of C, and C s , given 
by the relation above, 
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,(*> 



« 1 f 

(T^y 



(k + s)(k+s-1) „(* - 1) (*-<)(*-.*+ I) ^C* + 1) 



As 



As 



} 



sin 9 |3 



in which -77 -r^ or n _ tt »)» is ^4^ = 2,3015505 



37. Making * = -3, C + = g^ • ~x- . (C, -C 4 j. Using this 



formula, 

(*) 1 
C 3 =~rX 3,403041 

(5) 1 

C 3 = -r X 2,652559 

(6) 1 

C. = 7 X 2,047192 

(7) 1 

C, =-tX 1,568093 



(8) 1 

C, =-f X 1,193991 



(9) 1 

C 3 = -r X 0,904785 



(10) 1 

C,. = ^ X 0,682935 

(") 1 
C, = -3 X 0,513799 



(12) 1 

C, =-rX 0,385521 

T ft 

C 3 13) = K X 0,288655 



(14) 1 

C, = 7 X 0,215803 

(15) 1 

C; =-,X 0,161251 



(16) 1 

C, =~,X 0,120579 

(17) 1 

C, = -r X 0,090452 



„ rt ., , . 3 _,(« sin 9 f(2A + S)(2A+ 1) „(*-D 
38. Making. = -^, C^ = ^ . |* {£ C # 

_ v is c^ v . By the use of this formula, 

(S) 1 (9) 1 (13) 1 

C, = -r X 27,43922 C. =-rX 12,88246 C. =-,X 5,24565 

-s* ft ■§" u- -g- ft 



(6) 1 

C s =-tX23,14387 



(10) 1 (14) 1 

C; =-rX 10,39741 CI =-7X4,12790 

T ft y ft 



(7) 1 (11) 1 (15) 1 

C, = -r X 19,25046 C. =yX 8,32969 C. = ~/X 3,23120 

(8) 1 (12) 1 (16) 1 

C^ =^r X 15,82608 C f = 7 X 6,62955 C^ = ^ X 2,51561 

39 Makimr s - ± C W ^- / g*±*>M±i? C ( * ~ X) 

d9. Making * - fi , c { — 10cosa/3 • J ^ % 



(2&-5)(3&- 3) ^(*+i) 



C s }• . Thus we get 
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(6) 1 (9) 1 (12) 1 

C, =-rX 221,8780 C 7 =-, X 143,6296 C, = ~r X 84,9489 

^( 7 ) 1 ^( 10 ) 1 (13) 1 

C 7 =-rX 194,2735 C r =-, X 121,5988 C r = -r X 70,2184 

( 8 ) 1 (11) 1 (14) 1 

Cj. = ¥ X 167,9770 C, = Z X 102,0404 C, = ¥ X 57,6762 

(15) 1 

C } =^X 47,1003 
7 n (*) sin* /3 f (8j + 7) (2 & + 5) (* - 



40. Making . = I Cf = ^{ fl* ± "/" ± " C 

C 7 > . From this, 
C i' = ? X 1830 ' 596 C i 10) = 7 X 1266 ^09 C^ = yx 807,945 



(3£-7)(2£-5) ^(*+ 1) 



^( 8 ) 1 (11) 1 (14) 1 

C * -S X 1636,049 C + = Z X 1099,213 Cj = ? X 685,214 

( 9 ) 1 (12) 1 

C, =-rX 1446,655 C, =jX 946,016 

41 Making * - ^ C W - -^- J («* + 9)(8* + 7) c » - '> 
4i. waning * _ a , u v _ 18cos2|3 1 * . Vj- 

- j C, > . From this, 

( 8 ) 1 (10) 1 (12) 1 

C v = 7 X 15366,90 C v =^X 12473,68 C^ = -~X 9786,59 

(9) 1 (11) 1 (13) 1 

c v = Z x 139 °7 5 74 C v = ? xl 1092,76 c l v = 7 X 8570,07 

Section 13. 

Numerical calculation of (0,1) C s , (1,0) C s , 5fc. 

42. It will be sufficient to form, by differentiation, the expression for one of 

the differential coefficients of each order, as the others can then be derived by 

(£) 
simple addition. For C t is a function of a' and a of — I dimension : hence 

mdcccxxxii. o 
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a 5a' C * + ^ C * = " C s > or (^°) C s + i ' 1 ) C * = - C , • A g am ( as 

another instance) j-, 3 da C s is a function of a' and a of — 5 dimensions; con- 

d? W d 5 W d* U c ) 

seqnentlja' J ^ 7I - a C s +a^ KJ ^C s = -5^r^C s : or, multiplying 

(*) (*) (*) 

both sides by «' 3 a, (4,1) C; + (3,2) C; = - 5 (3,1) C; . It is indifferent 

which coefficient of each order we calculate first ; and for the algebraical pro- 
cess it is rather most convenient to differentiate successively with regard to 
the same quantity (as a'). 

ill i 



2 " "" (4 + |r-2co S% )' 

or, taking the coefficient of cos k% in the expansion on both sides, 

dor c « - a * a' u » + V " + «V * s+1 
Differentiating this formula with respect to a', and using the same formula to 

d? (*) 

simplify the differential coefficient, we get jj^ C s . In the same manner 
j-is C s , &c. are found ; multiplying them (beginning with j^, C s j by a, a" 1 , 
a' 5 , &c, we obtain the following expressions : 



<M» c?=- *«£> + (- 4- + .)..<f +1 

(2,0) c?'= +4 cf + (4 - »«) »•<?!,+ (- 4 + «) ! -^-JTT .c 

(M)cf— V 5 cf+-|(-T+»-)'-<f + . 



(*) 

s + 2 






J s + S 
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+ l(^-«)(-4+ 5a )- 5 -*~+" 1 -c 

/ 1 \3 (ft) 

+ \-^+V .«.* + i.* + 2.c; +8 

+ 2 (- 1 + „) 2 .(|-7«). *.r+i. 7+2. cfj 

/ l \4 (ft) 

+ \~^ + V -«.*+l.* + 2.* + 3.C f + 4 

(5,o ) cf' = -!icf + f(-i + 9 «) s .c« 1 

+ ?{-a-3«)(|-7«) + 4} S .F+l.C« 

+ t(t-»){-(t-3«)(|-7«)+4},.J+1.7+2.C 

K / 1 \3 / 1 \ (*) 

+ (--J + a ) . s . J+~T .* + 2.s + 3.* + 4.C 



(ft) 

s + 3 



-'s+4 



44. Using the same value of a as before, and making * = -% , these expres- 
sions become, 

(ft) l (*) (*) 

(1,0) Cj = - -g Cj - 0,3295790 . C, 

(2,0) C^ } = + | Cf - 0,3937533 . cf + 0,3258670 . cf 

(k) TS (ft) (ft) (*) (*) 

(3,0) Cj = - ^ Cj + 2,5134426 . C^ + 1,6567557 . C^ - 0,5369945 . C, 

o2 
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CO 105 CO CO O 

(4,0) C, = +-jq C t - 13,8031342 . <X - 6,6980500 . C^ 

CO CO 

- 5,9973139 . C r + 1,2388750 . C, 

CO 945 CO CO ^CO 

(5,0) C\ = - ~ C\ + 84,1230534 . C\ + 15,4872787 . C^. 

CO (*) „C0 

+ 10,2085636 . C, + 24,0925995 . C 9 - 3,6747654 . C, , 

T T "2" 

3 

45. Making s = -g, the formulee give 

CO 1 CO CO 

(1,0) C, = - -£ C; - 0,9887370 . C. 

CO 3 CO CO (*) 

(2,0) CI =+jC s - 1,1812599 . C, + 1,6293350 . C, 

CO 15 CO (*) (*) „C0 

(3,0) C, = - ~ C; + 7,5403278 . C. + 8,2837785 . C T - 3,7589615 . C 9 

5 

46. Making s = -g, the first formula gives 

V') 1 CO „V<) 

(1,0) c; = - ~ c; . - 1,6478950 . c; . 

(*) (*) 

47. Substituting in these the values of C. , C 3 , &c. found in the last sec- 
tion for different values of k, we form the following tables: 



For the development of the first term, 
k=8 (0,0)0^=^x0,0414571 

(l,0)C ( t 8) =7X -0,414243 
(2,0) cf =Ax4;i815 
(3,0) cf sr-W- 61,0595 
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(4,0) C ( x 8) = yx 897,236 

(5,0) C ( r 8) = Ax -14993,97 

* = 9 (0,0) C ( x 9) = — , x 0,0283925 

(1,0) C^ 9) = i X - 0,312394 (0,1) c|° = 1 x 0,284001 

(2,0) C (9) = -ix 3,S6300 (1,1) c£ 9) = ~ x - 3,23821 



3 <X 



i a' 



23361 



(3,0) C< 9) = ' x - 53,5643 (2,1) cf } = -L x 41,9753 

a a' i a' 

(4,0) C< 9) = lj x 832,244 (3,1) d 9) = -. X - 617,987 

i a' i a' 

(5,0) Cf } = J, x - 14512,93 (4,1) c[ 9) = ' x 10351,71 

i a' i a' 

k = 10 (0,0) C* 10) = i- X 0,0195495 

(1,0) Cf 0) = -I x - 0,234856 (0,1) c£ 0) = -I x 0,215306 

12 <Z 3£ ft 

(2,0) C ( r 10) =ix 3,13393 (1,1) C ( , 10) = X X - 2,66422 (0,2) C ( x 10) = 1x2, 

^ # so ^ a 

(3,0) Cj 0) = Jj x - 46,3921 (2,1) c£ 0) = -i X 36,9903 (1,2) cf 0) = J- x - 28,9976 

aft ^ fit -aft 

(4,0) d 10) = 1, X 761,088 (3,1) cf 0) = i x - 575,520 (2,2) C ( t 10) = ±, x 427,559 

(5,0) C ( ' 0) = -1 x - 13860,27 (4,1) cf 0) = JL x 10054,83 (3,2) C^ 0) = 2 x - 7177,23 
a a a a' A a' 

k = 11 (0,0) C ( r n) =-i x 0,0135189 

(1,0) C? 1} =ix- 0,176097 (0,1) C* U) = -l x 0,162578 

f a t a' 

(2,0) C^ U) =-lx 2,52220 (1,1) c£ U) = i x - 2,17001 (0,2) d n) =ix 1,84485 

s a' -s a' js a 

(3,0) 4 U) = -, x - 39,7288 (2,1) c[ ll) = -L x 32,1622 (1,2) C< U) = -I x - 25,6522 

a O a ' # a rt 

(0,3) C ( , 1I} = -i x 20,1 1 76 (4,0) cj. n) =-i X 687,024 . (3,1) C< U) = -L x - 528, 109 

^ fl a" CI "3 ft 

(2,2) d !1) = -1 X 399,460 (1,3) C< X1) = JL x - 296,851 (5,0) d U) = 1 x -13066,87 

E ft a" ft a ft 
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(4.1) c£ U> = JL x 9631,75 (3,2) cj- n) = -, X - 6991,20 (2,3) C* 10 = 1 x 4993,90 

k = 12 (0,0) C^ 12) = Ix 0,0093812 

(I,0)C ( ' 2) = Ix -0,131750 (0,1) C ( . I8) = 4 x 0,122369 

i a' i a' ' 

(2,0) C* 1S) = -I x 2,01559 (1,1) C^ 12) = -i- x - 1,75209 (0,2) C< 12) = 1 x 1,50735 

(3,0) C C _ 1S) = ~ x - 33,6822 (2,1) C< 12) = -^ X 27,6354 (1,2) C { * 2) = ' x - 22,3791 

(0,3) C ( x 12) = -L x 17,8570 (4,0) C<; 12) = JL x 612,866 (3,1) cf 2) = -L x - 478,137 

C£ Or 2 CI 7% (X 

(2.2) C ( x 12) = -1 x 367,595 (1,3) C< 12) = -L x - 278,079 (0,4) C { ™ } = _L x 206,651 

(5,0)C ( ' 2) =-Lx -12169,39 (4,1) C< 12) = -1 x 9105,06 (3,2) d 12) = 1 x -6714,37 
i a' 2 a' i a 

(2,3) C ( T 12) = -1 x 4876,40 (1,4) C ( x 12) = -1 x - 3486,00 

A = 13 (0,0) C ( T 13) = -i- x 0,0065274 

5 a 

(1,0) C^ 13) = -^ x - 0,098398 (0,1) C< 18) = i- x 0,091871 

2 CI a # 

(2,0) C ( x 13) = -1 x 1,60062 (1,1) 4 1S) =ix- 1,40382 (0,2) C< 13) = -I x 1,22008 

7£ Ct 2 CI % Q. 

(3,0) C ( T 13) =Ix- 28,3028 (2,1) C ( z 13) = ~ x 23,5009 (1,2) C< 13) = ~ x - 19,2894 

(0,3) C ( r 13) = -I x 15,6292 (4,0) C^ 3) = -l x 540,744 (3,1) C^ 3) = -1 x - 427,533 

i a' £ a # «' 

(2.2) C ( r 13) =1- x 333,529 (1,3) cf = -L x - 256,371 (0,4) C^ 3) =-L x 193,854 

(5,0) C^^-l x - 11205,35 (4,l)d 13) =-L x 8501,63 (3,2) C< 13) = JL x -6363,96 

i a £ a £ a' 

(2.3) C^ 13) = J- x 4696,32 (1,4) C^ 13) = 1 x - 3414,46 (0,5) C^ 13) = -I x 2445,19 

For the development of the second term, 

k = 9 (0,0) C (9) = -1 x 0,904785 

•|- a 

(l,0)C (9) =-ix -13,18976 
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(2,0) C (9) = -1x219,4819 
4- a 

(3,0) C (9) = -L x -4152,686 

-§■ a' 

k = 10 (0,0) C (10) = -1 X 0,682935 

(1,0) C (10) = 1 x - 10,62177 (0,1) C (10) = 1 x 9,93883 

(2,0)C (1O) =-L x 186,3554 (1,1) C (10) = -1 x -165,1119 

(3,0) C (I0) = i, x - 3677,095 (2,1) C (10) = -1 x 3118,029 

4 «' -f « 

4=11 (0,0) C (U) = 1x0,5 13799 

■f a' 

(1,0) C (H) = 1 x - 8,49277 (0,1) C ( , H) = ± x 7,97897 

(2,0)C (,1) = i x 156,8038 (1,1) C (H> = 1 x - 139,8183 (0,2) C (U) = -1 X 123,8604 

-§• a -§• a -f a 

(3,0) C (U) = -L x -3224,776 (2,1) C (11) =-1 x 2754,365 (1,2) C (U) = -1 x -2334,910 

■f a, -f a t a 

A = 1 2 (0,0) C (12) = -L x 0,385521 

-§• a 

(1,0) C (12) = -L x - 6,74764 (0,1) C ( 3 12) = i X 6,36212 

(2,0) C (12) = 1 x 130,8682 (1,1) C (12) = -1 x - 1 17,3729 (0,2) C ( / 2) = ~ x 104,6487 

-§■ (X -^- fit -j- w 

(3,0) C° 2) = i X - 2803,076 (2,1) C (12) =-L x 2410,471 (1,2) C^ = ~ x -2058,352 
j. a ■§■ & t a 

(0,3) C (12) = ~ x 1744,406 
■J- o 

For ^e development of the third term, 

k = 10 (0,0) C (10) six 10,39741 

(l,0)C ( ' O) = ix -205,5808 

A =11 (0,0) C (U) = -1x8,32969 

(1,0)C ( 5 U) = -1 x -172,3167 (0,1) Cj" 5 ^ x 163,9870 
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(8) (9) (9) 

48. Employing these numbers in the calculation of L , L , &c. M , 

(10) (10) (11) 

M ,&c. N and N , from the expressions in (20), (24), and (28), we 
obtain the following numerical values : 

( 8 ) m 

L = —, X - 333,0969 

(9) m 

L =-^X 1273,4929 

(io) m 
L =-, X -1945,7913 

(ii) m 
L =-,X 1485,3152 

T (is) m 

L =^X -566,5632 

(13) m 
L =5>< 86 ' 3635 

(9) m 
M = —, X - 503,4795 

(io) m 
M =-,X 1088,9148 

(H) m 
M K = 5 X - 787,0581 

(12) m 
M K =|x 190,0487 

(io) m 
K =^X -85,3347 

(ii) m 
N = -, X 58,8603 

49. The computation of these quantities has been effected by means of alge- 
braical operations of great complexity, and numerical calculations of no in- 
considerable length ; and it is not easy to find in the operations themselves 
any verification of their accuracy. This has imposed on me the necessity of 
examining closely every line of figures before I proceeded to another. I have 



IN THE MOTIONS OF THE EARTH AND VENUS. 105 

had the advantage however of comparing the calculated values several times 

with the values which I calculated nearly four years ago. At that time I 

developed the principal fraction in a different manner, and I expressed the 

CO 
quantities C„ &c. by different formulae ; and the fundamental number differed 

by a few units in the last place of decimals. The numbers admitted of com- 
parison at several intermediate points before arriving at the final results ; and 
one small error was discovered in the old calculations, and one in the new 
ones. Upon the whole, I am certain that there is no error of importance in 
these numbers ; and I think it highly probable that there is no error, except 
such as inevitably arise from the rejection of figures beyond a certain place 
of decimals. It is impossible to assert that the last figure preserved is correct, 
or even the last but one ; but 1 do not think that the last but two is wrong. 



Section 14. 

Numerical calculation of the long inequality in the epoch, depending on 
(13 X mean long. Earth — 8 X mean long. Venus). 

50. The most convenient form in which the expression of (29) can be put is 
the following. 

| L (8 . e' 5 . cos (5 vf) + L (9) . e* e . cos (4 ti + w) + L (10) . <?' 3 e 2 . cos (3^+2 m) 

C 11 ) , . (12) 

+ L .*' 2 « 3 .COS(2w' + 3w) + L . e' e 4 . COS (sr' + 4 sr) 

(13) (9) 

+ L . e 5 . cos (5 a) + M . e' 3 / 2 • cos (3 w> + 2 6) 

+ M (10) . e' 2 ef 2 . cos (2«r' + sr + 2 6) + M (U) . e' e 2 / 2 . cos (m> + 2 v + 2 S) 

(12) (10) 

+ M .e 3 / 2 .cos(3sr + 20)+N . e'/* . cos (a* + 4 6) 
+ N (U . e/ 4 . cos (w + 4 &) j cos {13 (n't + s 1 ) — 8 (n t + s) } 

+ | L (8> . e* . sin (5 *j>) + l} 9 \ ^ e . sin (4 W + w) + L° 0) . e s e 2 . sin (3^+2 ») 

(ii) ri2) 

+ L .e' 2 e 3 .sin(2w' + 3w)+L . e' e i . sin (w> + 4 w) 

MDCCCXXXII. P 
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+ L (13) . £ . sin (5 «) + M (9) . J*j* . sin (3 W + 2 6) 

+ M (10) .^ 2 e/ 2 .sin(2«7'+ w 4-2^)+M (11) ./e 2 / 2 .sin( C r'+2tif+2^) 

+ M (12) . e3/2 . sin (3 « + 2 0) + N (10) . </» . sin (^ + 4 6) 

+ N (U) . <?/ 4 . sin (ar + 4 0) j . sin {13 (n'* + s') - 8 (»* + «)} 

The elements <?', <?, &c. are all subject to small permanent variation ; and 
(considering the great length of period of the inequality which we are calcu- 
lating,) those variations may have a sensible influence upon it. It is prudent 
therefore; as well as interesting, to take into account these variations. 

5 1 . Let P and Q be the values of the coefficients of cos {13 (n't + s') — 8 (nt + 1) } 
and sin {13 (n't + «') —8(nt-\- s)} in the expression above, giving to the ele- 
ments the values which they had in 1750. Then, as all the permanent varia- 
tions are small, the powers of t above the first may be rejected, and the coeffi- 
cients at the time t after 1750 may be represented by P +pt andQ + q t. 
Thus the term of R becomes 

(P+pt) cos {13 (rit + i') —8 (nt + s)} +(Q + q t) sm {13 (n't + i 1 ) — 8 (nt + s)} ; 
and by (2), omitting the terms depending on j^ and -^p for the reasons in (31), 

g=-^f?(P +| ,*) B ln{l3(«f# + 0-8(n* + i)} 

+ ^7^-' (Q + ? *) cos { 13 (n't + - 8 (n t + »)} 



dt ~ ^ <*' 

39 n' s a' 



(Qt + q t 2 ) cos { 13 (n' t + - 8 (n t + t) } 



Integrating these, (considering »', s', n, and e, on the right-hand side, as 
constants,) and substituting in the expression n't + g', it becomes 

, 39ri*a' f P +p.t Qq 1 , , ,. a / , , n, 
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The terms added to N' t + E' constitute the inequality in the epoch. 

52. The values of the elements for 1750 and their annual variations are 
given by Laplace in the M6canique Celeste, 2 me Partie, Livre 6, N os 22 and 
26. To give them the form necessary for our purpose, we must from the varia- 
tion in a Julian year deduce the variation for a unit of time. Now a Julian 
year is (nearly) the time in which the angle ri t increases by 2 it ; its expression 

is therefore —r. Consequently if we multiply the annual variations by g— , we 

shall have the variations in a unit of time : and if we multiply them by „— , we 

shall have the variations in the time t. With regard to the quantities \j}, &c. 
introduced by Laplace for the purpose of altering his assumed masses if neces- 
sary, it may be observed that the only planet which materially affects the 
changes of the elements, and whose mass is known with certainty to require a 
change, is Venus herself. The investigations of Burckhardt and Bessel lead 
to the same conclusion as my own (Phil. Trans. 1828), namely, that the mass 

8 1 

of Venus is -5- X the mass assumed by Delambre, or 401211 of the sun's mass. 
Laplace supposed it 3 B3137 of the sun's mass : the comparison of these gives 

Laplace's \j} = — ,045. In using Laplace's expressions, therefore, I shall sup- 
pose (a' = — ,045, and /*, p", \j}'\ &c = 0. For convenience, the centesimal * 
division will be retained. 

53. Thus we have 

_ 650198000 . , 
71 ~ 399993009 X U 

e' — 0,01681395 — 0,0000000729 X fit 

e = 0,00688405 — 0,0000001005 f X n't 

f = 0,02960597 + 0,0000000172 X n't 

s/= 1098,5790 + 0,0000091017 X n't 

* Bokda's tables, published by Delambre, have been used in these computations. 

f The variations of the elements of Venus do not agree with those of Lindenau's tables. 

p2 
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sr = 142&,1241 — 0,0000018080 X lit 
6 — 82S,7093 — 0,0000139997 X ri t 

The node and inclination are those on the earth's true orbit. All the coefficients 
of ri t are in decimal parts of the radius 1, and not in parts of a degree. 

54. From these we deduce the following values, the figures within the 
brackets being the logarithms of the numbers. 

e' 5 = + (91,1283485) — (86,46438) . ri t 
e' 4 e = + (90,7405229) — (86,24488) . ri t 
e'3 e 2 = _j_ (90,3526973) - (85,97806) . ri t 
e'2 e 3 = + (89,9648717) — (85,68477) . ri t 
e'(A — + (89,5770461) - (85,37453) . ri t 
e* = + (89,1892205) — (85,05252) . ri t 
e '3f2 _ + (91,6197677) - (86,69331) .ri t 
e' 2 ef 2 = + (91,2319421) — (86,57650) . ri t 
e ' e 2/2 — -j- (90,8441165) - (86,35426) . ri t 
e 3 / 2 = + (90,4562909) - (86,08606) . ri t 
<//4 = + (92, 1111 869) - (86,4 1 479) . ri t 
ef* — + (91,7233613) — (86,81239) . ri t 

( 8 ) 

a ' L . cos (5 «/) = + (2,3572098) + (98,04404) . ri t 

771 
(8) 

t-1— . sin (5 */) = - (2,3859510) + (98,01530) . ri t 

Tit 

rr (9) 
a lj cos (W + ») = - (3,0841670) - (98,12469) . rit 



m 

(9) 



d L . sin (4 t*' + w) = + (2,5856285) - (98,62323) . ri t 



m 



, , (io) 
a tJ cos (3 ra-' + 2 w) = + (3,2799989) - (97,97020) . ri t 



a'L 



m 

(10) 



m 



. sin (3w' + 2w)= + (2,5956493) + (98,65455) . ri t 



— — . cos (2 ra' + 3 w) = - (3,0497482) + (98,09507) . ri t 
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I T (11) 

— sin (2 vf + 3 m) = - (2,9885633) - (98,15626) . n t 

Tib 

(12) 

a !k_ . cos («/ + 4.w) = + (2,2816808) - (96,99874) .n't 

(12) 

a ' L . sin (W + 4 w) = + (2,7269678) + (96,55345) . ri t 

m 

, T (13) 

"l}±— . cos (5 sr) = + (1,1565787) - (96,88643) . ri * 



. sin (5 sr) = - (1,9302586) — (96,11275) .m'^ 



/ T (1S) 

to 

(9) 

a ' M cos (3sr' + 2^) = - (1,6642318) — (96,54170) . ri t 



m 



a'M 



a' 



in 
'M (10) 



. sin (3 vt' + 2 0) = - (2,7001492) + (95,50578) . »' t 



m 



cos (2 w' + w + 2 0) = - (2,6468283) + (98,06223) . ri t 



,(i°) 
^M__ . sin (2 <*' + w + 2 0) = + (2,9976206) + (97,71 143) . ri t 

7ft 

(ii) 

— cos (V + 2 *7 + 2 0) = + (2,8001796) - (98,02467) . ri t 

(ii) 
— — . sin (•» + 2 w + 2 4) = - (2,6722198) - (98,15263) . rc'* 

TO 
(12) 

— cos (3 vr + 2 ff) = - (2,2752440) + (96,91212) . «' * 

7/6 

(12) 

— — . sin (3 «r + 2 0) = + (1,3880761) + (97,79929) . ri t 

(10) 

fflfN . COS (ot' + 4 6) — — (1,8370062) - (97,37538) . ri t 

Tilt 

fM (10) 

— sin (vs* + 4 6) = - (1,7042256) + (97,50816) . »' * 

77i 
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— — . cos (w + 4 6) = + (1,3847917) + (97,49139) . ri t 

7fl 

— sin («■ + 4 0) = + (1,7294138) - (97,14677) .n'* 

55. Substituting these in the expressions of (50), we find 

P = _ -J X (94,1302623) p = + ~j X (89,08397) . ri 

Q = - -J x (94,0722348) g = + -^ X (89,47976) . »' 

j i • ? " l j !o r o 1674883 v> , . 

and making -7- = 4Qian , and 13 w' — 8 n = — 309903090 X », in the expression 

of (51), we find for the long inequality 

{ — (94,8787039) + ri t X (89,82780)} . sin { 13 (ri t + g') - 8 (» £ + s) } 
+ {+ (94,8139258) — ritX (90,22359)} .cos {13 (»'* + s') — 8 (n* + 1)} 

which may be put in the form 

{ + (94,9992364) - ri t X (90,20461)} . sin {8 (nt + s) - 13 (ri i + g') 

+ 40° 44' 34" — ritX (94,91918)} 

where the degrees, &c. in the argument are sexagesimal. The coefficient is ex- 
pressed by a multiple of the radius : to express the principal term in sexage- 
simal seconds, it must be divided by sin l". And if Y be the number of years 
after 1750, since rit — mean motion of the earth in Y years = 2 t . Y~6 .60 3 . Y 
in seconds, the coefficients of rit must be multiplied by 6 . 60 3 . Y, and their 
values will then be exhibited in sexagesimal seconds. Thus we find at length 
for the inequality 

{2",059 - Y X 0",0002076} X sin {8 (n t + 1) - 13 (ri t+s') 
+ 40° 44' 34" - Y X 10",76}. 

56. The mean longitudes n t + &, ri t + s', are measured from the equinox of 
1750. But if I, V, are the mean longitudes of Venus and the Earth measured 
from the place of the equinox Y years after 1750, then (in consequence of pre- 
cession) 
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nt + s=zl—Y X 50",1 
lit + s = I'- Y X 50", 1 

Consequently 8 (n t + s) — 13 (n't + s') = 8 / - 13 /' + Y X 250",5. 

Substituting this, the expression for the inequality is 
{2",059 - Y X 0",0002076} X sin {8 I - 13 I' + 40° 44' 34" + Y X 239",7> 

57. I have compared the calculations of the principal part of this inequality 
with the calculations made in 1827. Two errors were discovered in the former 
calculations, one of which was important. I am quite confident that there is 
no sensible error in the results now presented. The terms depending on Y 
were not calculated on the former occasion : but the calculations now made 
have been carefully revised. 

Section 15. 
Numerical calculation of the long inequality in the length of the axis major. 

58. This being very small, we shall omit the variable terms. Thus we have 

~- Q . cos {13 (n't+ s')~8(nt + $)} 



P- 



whence 



= A' - a' . (92,31993) . cos {13 («' t + — 8 (» t + a)} 
- a' . (92,26190) sin {13 (ii t + s') - 8 (n t + s) } 

= A' - a' X 0,000000027756 X cos {8 (« ^ + 1) - 13 (%' t + e') + 41° 1 1'} 

The magnitude of the coefficient is barely ^Vth of Laplace's minimum, and 
this inequality may therefore be neglected. 
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Section 16. 
Numerical calculation of the long inequality in the longitude of perihelion. 

59. The expression for -jj being rg . jj, the part which we have to 

consider may be put under the form 

dts' rid f cT ( 8 ) ,- ,_ , N . , T (9) ,. ,. , , . 

-jj- = — ■ -7-^5 . -j 5 L . e' 5 . cos (5 ro-') + 4 L . e 4 <? . cos (4ot + cr) 

(i°) . (") 

+ 3 L . e' 3 e 2 . cos (3 «/ + 2 w) + 2 L . e 2 e 3 . cos (2w'+3 «r) 

(12) (9) 

+ K . 4* . cos («/ + 4 w) + 3 M . e' 3 / 2 . cos (3 m + 2 4) 

(10) (11) 

+2 M . <?' 2 e/ 2 . cos (2 w + w + 2^ + M . <?' e 2 / 2 . cos (w'+2 *r+2 4) 
+ N (10) . e'/ 4 . cos (w + 4 | cos 1 13 (»' * + /) - 8 (n * + s) } 

M'ffi' f ( 8 ) (9) 

~" T7* i 5 ^ • e ' 5 • sm (5 ^') + 4 L . e' 4 e . sin (4 w' + vr) 

(10) (in 

+ 3L ./ 3 e 2 .sin(3s/+2sr) + 2L . e' 2 e 3 . sin (2 *' + 3 «r) 

(12) (9) 

+ L . e' e 4 . sin («•' + 4 «■) + 3 M . e' 3 / 2 sin (3 «/ + 2 6) 

+ 2 M (10) . e' 2 ep . sin(2w'+w + 2 ^ +M (H) .^ 2 /2. sin(w'+2w + 20 

+ N (10) . e'/ 4 . sin (70-'+ 4 4) jsinj 13 (»'* + g') - 8 (nt + «) j 

which (neglecting the variable terms) is found to equal 

n X (92,35866) . cos { 13 (n' t + e') - 8 (» * + e) } 

+ »' X (92,60190) . sin { 13 (n' * +• «') - 8 (» t + «)} 

Integrating, 

w = n' - (94,73673) . sin { 13 (n * + s') - 8 (n * + s) } 

+ (94,97997) . cos {13 (»' * + - 8 (n * + s)} 

or 

w = IT + 1",1250 . sin {8(nt + s}- 13 {tt t +g')} 
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+ l",9697 . cos {8 (n t + s) - 13 («' * + s')} 
= IT + 2" 2683 . sin [8 (n t + i) - 13 (n t + s') + 60° 16'} 

Section 17. 
Numerical calculation of the long inequality in the excentricity. 

60. On forming the expression for -j^, or + -zrji • j^,, it is immediately seen 

that the coefficients of cos {13 (n't + /) — 8 (nt + s)} and sin {13 (n't -f- s') 
— 8 (nt + s)} are related to those above, and that 

^j = + eV X (92,60190) .cos {13 (V £ + e') — 8 (nt + «)} 

— erix (92,35866) . sin { 13 (n t + e') — 8 (n t + s)} 

Integrating, 

e' = E' - e X (94,97997) . sin { 13 (it t + t) — 8 (n t + 1) } 

- <?' X (94,73673) . cos { 13 (n' t + e') — 8 (nt -{- s)} 

= E' — (92,96240) cos {8 (n t + s) — 13 («' * + g')} 
+ (93,20564) . sin {8 (n t + s) — 13 (rc'tf + $0} 

= E' - 0,0000001849 . cos {8 (n t + s) — 13 («' * + a') + 60° 16'} 

The principal inequality in the radius vector is that produced by the last 
term : it is however too small to be sensible. 



PART II. 

PERTURBATION OF THE EARTH IN LATITUDE. 

Section 18. 
Explanation of the method used here. 

61. If <f>' be the inclination of the earth's orbit to the plane of xy, and (f the 
longitude of the node, then 

mdcccxxxii. Q 
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il — n ' a ' L ^5 

dt~ ~ pf 4/T^Y* ' <P' ' df' 

d ^ . n ' a ' _L ^5 

dt — ~T~ pi t/T^Z 2 ' i ' dtf 

or, neglecting e' 2 , 

dV _ n'a' J_ dR 
dt~ p! ' ' <p' ' d <ft 

d<~" + a' m <ff ' dV 

These expressions are true only when <p' is so small that its square may be re- 
jected. This restriction, however, is convenient as well as necessary. For in 
the expansion of R we shall have to proceed only to the first power of <p', and 
make <p' = when we have arrived at our ultimate result : consequently the 
same values of 6 and <p must be employed as in the first Part. 

62. The only term of R, which by expansion will produce terms of the 
form cos (13 — 8) with coefficients of the fifth order, is the fraction 

— — TYh 

' */{(!*! "— x ? + (iS—vY + (z' — zJ H' ^ or SUDS titution in the denominator we have 

x = r cos v' (neglecting <p 2 ) 

y z= r sin v 

z = r . <p' . sin (v — &') 

x = r {cos (v — ff) . cos 6 — cos <p . sin (v — 6) . sin 6} 

y — r {cos (v — 6) . sin 6 -\- cos <p . sin (v — 6) . cos 6} 

z =r . sin <p . sin (v — 6) 

whence the fraction is changed to 

— m 

\/ {i n — Qr'r. cos(rf -v)+r i + Qi J r.f*. cos(rf — v)— 2rV./ a .cos(z/ + w— 2fl) — 4r / r.<p[/.sin(i/— S').sin(u— 0)} 

where/ is put for sin -^ and 2/ for sin <p, on the principle of (13). The part 

of this depending on the first power of <p' is 

—m.Qr'r. <p'f. sin (p* — fl') . sin (v — 6) 

{r' 2 — 2 r' r . cos (p' -») + r s + 2r l j' .f" 2 . cos (p' — v) — Ir 1 r ./* . cos (p' + v — 2 0} T 
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of which, on the principle of (8), &c, we are to take only 

m .r' r. <p'f. cos (t/ + v — 6' — fl) 
{r 12 — 2 r 1 r . cos (w' — v) + r* — 2 r 1 r .f* . cos (»' + w — 2 0)} 4 

Expanding the denominator by powers of/" 2 , this becomes 

m.r'r.tff. cos(i/+ p — $' — 6) w . 3 r 12 r 3 . $y* . cos (■d + v — V — fl).cos (p' + P— 2fl) 
{r' s — 2 r' r . cos (t/ — v) + r 8 }* ' {r' 2 — 2 r'r. cos (p' — p) + r a }* 

or 

ot . ?•' r . §'f. cos (p' + g — fl' - ^) , 3 ot . ?-' 3 r a . fl'/ 3 . cos (2 p' + 2 p — $' — 3 fl) 
|,/3 _ 2 / r . cos (p' — v) + r 3 }^ 2 {r /2 — 2 r r . cos (xf — p) + r 9 } 1 " 

Section 19, 

Selection of the coefficients of cos (13 — 8) in the development of the two last 

fractions. 

63. If we compare the first fraction with the fraction developed in Section 8, 
we perceive that the following are the only differences between them. The 
signs of the coefficients are different : and in the coefficient of the new fraction 
(and in every term of its development) there is <p' instead oif with the corre- 
sponding change of argument. From this it is readily seen that the coefficient 
of cos (13 — 8) will be formed from that in Section 8 (Art. 24), by changing the 

a' 
sign and multiplying by y ; the argument always being changed according to 

the rules of (9). The coefficient is therefore 

- M (9) . e' 3 <p'f - M (10) . e' 2 e <p'f - M ( "\ e' e* <p' f - M° 2) . e 3 <flf. 

64. If we compare the second fraction with the fraction developed in Sec- 
tion 10, we see that there are the same differences as those mentioned above, 
with this additional one, that the multiplier is double of the multiplier of the 
fraction in Section 10. Thus the coefficient of cos (13 — 8) is found to be 

- 2 N (10) . e' <p'/ 3 - 2 N (n) . e <p'f* 

The sum of the terms in these two sets, multiplied respectively by the cosines 
of their proper arguments, constitutes the whole term of R which we have to 
consider. 

q 2 
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Section 20. 
Numerical calculation of the perturbation in latitude. 

(9) 

65. The first of the terms found in the last section is — M . e' 3 <p ! f. cos 
{ 13 (n' t + i) — 8 in t + s) — 3 -d — § — 6}. With respect to this term only, 

j? = - M (9) . e' 3 /- cos { 13 (n' t + a') - 8 (n t + s) - 3 «' - 4'- 0} ; whence 

(9) 

— 8 («* + f) - 3 «' — tt — 0}. 
And 

^§ - - M (9) . e 3 </. sin {13 (»' * + i) - 8 (»< + i) - 3 vt - 6'- 6} ; 
whence 

* =* + j??./ rff = * + -7- • 13n'-8n e3 /-COS {18(rf* + 

— 8 (» if + s) - 3 vt — $ - 6}. 

The Earth's latitude, neglecting small terms, is <p' . sin (n' t -f- 1 — tf). And 
from the expression above, sin (»' /+$' — /) = 

(9) 

M a' w' e' s f 

sin (rf, + ,'_ 0') _— j- . v3n ,_ 8n . -f.co&(ftt+<- 0').sin 

{13 (rc'if + — 8 (re* + «) - 3 <*' — • — 6} 

Multiplying this by the expression for <p', and putting /, <p\ for 0', <3>', in the 
small terms, we find for the latitude 

O^ s in(^^^ + ^ '-0O-^. I 3^^ i .e' 3 /•sin{12( ? ^^ + O 

— 8 (n t -f g) - 3 w — 6} 

and the last part, or the perturbation in latitude, is 

(9) 

n' M a' 

~ isF^sli ' e ' 3 /- -7- -sin {12 (»'* + ~ 8 (» * + f) - 3 «r' - 4} 

Similar expressions will be obtained from all the other terms. 
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66. If we put for sin { 12 (re' t -f- s') — 8 (re t -f- s) — 3 ■&' — 0} its equivalent 
cos (3 «/ + 2 0) . sin {12 (re' # + {) — 8 (re * + e) + 0} - sin (3™'+20). cos 
{12 (re'/ -f- «') — 8 (nt -f- g) + 0}, and similarly for the other terms, we find for 
the whole coefficient of sin { 12 (nt -j- %) — 8 (re t + g) + 0)}, 

— T5~7 5~ .— .-?. — "S e / • M . COS (3 tsr + 2 0) 

+ e 2 e/ 2 . M (10) . cos (2 n/ + « + 2 6) + e' e 2 / 2 . M (U) . cos (s/ + 2 « + 2 0) 

(12) (10) 

+ e3/ 2 .M .cos(3tf + 2 0) + 2<?'/ 4 .N .cos(«r' + 40) 

+ 2 <?/ 4 . N (U) . cos (w + 4 4) j 
and for the whole coefficient of cos { 12 (re' t + s') — 8 (re * + e) + 0}, 
+ T5F=8i- ?-7 ^{^/ 2 .M (9) .sin(3.'+2^) 

+ e' 2 e/ 2 . M (10) . sin (2»' + tf + 2H«' e 2 / 2 . M (U) . sin ( i /+2»+20) 
+ e 3 / 2 . M (12) . sin (3 v + 2 4) + 2 e'/ 4 . N (10) . sin (</ + 4 0) 



+ 2e/ 4 .N (U) .sin(sr + 4tf)l 



On performing the calculations, the inequality is found to be 
+ 0",0086 . sin {8 (» / + s) — 12 (re' t + g') - 0} 

+ 0",0060 . cos {8 (re * + s) - 12 (re' / + — 6} 
or + 0",0105 . sin {8 (re * + g) - 12 (re' * + g') - 39° 29'} 

which is too small to be sensible in any observations. 



PART III. 

PERTURBATIONS OF VENUS DEPENDING ON THE SAME ARGUMENTS. 

67. If we consider Venus as disturbed by the Earth, and take the orbit of 
Venus for the plane xy, the term involving cos (13 — 8) in the expansion of R 
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will be exactly the same as when we consider the Earth disturbed by Venus. 
For the longitudes of perihelia and the longitude of the node will be the same : 
the sign of f will be different, but as the even powers only of this quantity enter 
into the expansion of R, and as its magnitude (without respect to its sign) is 
the same in both, that circumstance makes no difference. It is only necessary 
then to put m' instead of m in the multiplier of the term. 

68. First, then, for the inequality in the epoch. Observing that in the ex- 
pression of (51) the multiplier m is included in P, p, Q, and g, it will be seen 

that for the perturbation of Venus we must use the multiplier in- 

stead of \ — -. That is, the argument of the perturbation of Venus is the 

same as that of the Earth ; and its coefficient is found by multiplying the co- 
efficient for the Earth by — 7ir , 2 * . As 8 n = 13 n' nearly, this fraction 

1 <J ? I Oj Jib " 

nam! 13 a m' , . m! 1 , , 

= -¥-^-m = -"8'-^-m nearl y- Assuming - = ^9630* and the other 

13 401 21 1 

quantities as before, this fraction is — -g- X 0,72333 X 329630 : wnencetne l° n g 
inequality in the epoch of Venus = 

{— 2",946 + Y X 0",0002970} X sin {8 I - 13 I' +40° 44' 34" + Y X 239",7} 
= {2",946 — Y X 0",0002970} X sin {8 J — 13 /' + 220° 44' 34" + Y X 239",7} 

The corresponding inequality in the axis major, like that for the Earth, is in- 
sensible. 

69. For the long inequality in the longitude of perihelion. This cannot be 
deduced from that of the Earth: but, calculating it independently in the same 
manner, it is found that 

sr = n - 0",008 sin (8 / - 13 V) - 5",704 . COS (8 I - 13 ff 

70. For the long inequality in the excentricity. This may be derived from 
that in the longitude of perihelion in the same manner in which it was done 
for the Earth : thus it appears that 

e = E — 0,0000001904 . sin (8 1 — 13 I') + 0,0000000003 . cos (8/ — 13 /') 

71. For the inequality in latitude. The orbit of Venus must now be sup- 
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posed to be inclined at a small angle to the plane of x y. We have remarked 
that, in the development of R for Venus as the disturbed body, the sign of/ 
will be changed : and as the term of R on which the perturbation in latitude 
depends is a multiple of odd powers of/, the sign for Venus will be different 
from that for the Earth, Besides this there will be no difference, except that 
a m' n is to be substituted for a' m n'. Proceeding then as in (65), and con- 
sidering the effect of the first term of (63), we find 

6 = — rr-j — — . — . — . . V —J- . sin {13 (n't + g') 

— 8 (» t + s) — 3 tf - 6' — 6} 
whence 

sm (nt + i - 6) = sin {nt + s - 0) + i3w >_ 8w • ~ • j • -£~ • / X 

cos (rc* + s — 0) .sin {13 (»'/ + 0-8(»« + 5)-3s'-^-^} 
And 

(9) 

p= O - ttt-j 5- . — . — . — — - . e' 3 /-COS{13 (n * + s) 

- 8 (» if + s) — 3 w — 4' — 6} 
The product of these expressions gives for the latitude of Venus 



m 



a M (9 V 



— 9 (nt-\-s) —3*x' — d'} 
where & has the same value which had in the investigation for the Earth. 

(9) 

The perturbation in latitude is therefore jj^rzTs^ . — ..— . e' 3 / . sin 

{13 (rc'tf + &') — 9 (rc* + s) — 3w' — ^}, and similarly for the other terms. 
Comparing this with the term in (65) it will readily be seen that we have only to 

yn M ft 

multiply the expression of (66) by — - —, —„ and to put 9 (rctf + g) — 13(»'*-f-s') 

instead of 8 (n t + s) — 12 (ri t + s'), and the perturbation of Venus in latitude 
will be found. Thus it becomes 
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— 0",0123 . sin {9 (n t + g) — 13 {n't + g') — 4} — 0",0086 . cos {9 (» * + «) 

- 13 (»'/ + - ^} 
or 

+ 0",0151 X sin {9(ra* + g) — 13(«'* + g') + 140° 31'} 

which, though larger than the Earth's perturbation in latitude, is too small to 
be observable. 

Conclusion. 

It appears, then, that in calculating the Earth's longitude (or 180° + Sun's 
longitude), the following terms should be used in addition to those that have 
hitherto been applied ; (where / and /' are the mean tropical longitudes of Venus 
and the Earth, and Y the number of years after 1750 :) 
To the epoch of mean longitude 

+ {2",059 - Y X 0",0002076} X sin {8 1- 13 V + 40° 44' 34" + Y X 239",7} 

To the epoch of longitude of perihelion 

+ 2",268 X sin {8 1- 1 3 I' + 60° 16'} 

To the excentricity 

- 0,000000 1 849 . cos { 8 1 - 1 3 /' + 60° 1 6'} 

and that, in calculating the Earth's latitude (or the Sun's latitude with sign 
changed), the following term should be used ; 

+ 0",0105 . sin {8 1— 12 t - 39° 29'} 

Similarly, it appears that in calculating the place of Venus, the following terms 
should be applied : 

To the epoch of mean longitude 

+ {2",946— Y X 0",0002970} X sin {8J- 13^ + 220° 44' 34" + YX 239",7} 

To the longitude of perihelion 
-5",70.cos{8/-13/'} 

To the excenti-icity 

- 0,000000190 . sin {8 /- 13 ?} 

To the latitude 

+ 0" 3 0151 .sin {9^-13/'+ 140° 31'} 
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The terms affecting the latitude may be at once neglected. The inequali- 
ties in longitude produced by the change of mean anomaly and excentricity, 
(n t + i — 73 and e), and which are 

for the Earth 

— 0^0470 X sin {8 1- 12 f- 15° 34'} -0",0346 . sin{14/'- 8 1- 139° 22'} 
for Venus 

+ 0" 3 0671 . sin{9 1- 13 /' — 24° 40'} +0",0203 . sin{13 I'-jl- 168° 40'} 

can scarcely be detected from observation. The inequalities in the radii vec- 
tores are not sensible. 

The long inequalities in the epoch of longitude are however by no means 
to be neglected. To point out a single instance in which their importance will 
be sensible, I will estimate roughly their effect on the places of the Earth and 
Venus at the next transit of Venus over the Sun's disk (in 1874). The value 
of these inequalities at the time of Bradley's observations was small ; and they 
were at their maximum at the beginning of this century. If, then, the mean 
motions of the Earth and Venus were determined by comparing the observa- 
tions about Bradley's time with the observations a few years ago ; the Earth's 
longitude in 1874, when the inequalities are nearly vanishing, would be too 
small by nearly 4" ; that of Venus would be too great by 6" : their difference 
of longitude would therefore be nearly 10" in error ; and this would produce 
on the geocentric* longitude of Venus an effect of between 20" and 30". As 
another instance, I may mention that the secular motions of the Earth, deter- 
mined from observations of two consecutive centuries, would differ nearly 8", 
and those of Venus nearly 12". 

These inequalities vanish in the years 1622, 1742, and 1861 ; and have their 
greatest values, positive for the Earth and negative for Venus, in 1682 ; and 
negative for the Earth and positive for Venus, in 1802. At the principal 
transits of Venus their values are as follows : 

* In the Memoirs of the Astronomical Society I have pointed out the utility of observations of 
Venus near inferior conjunction for determining the coefficient of the inequality in the Earth's motion, 
produced by the Moon. I take this opportunity of repeating my conviction, that observations of Venus 
near inferior conjunction are adapted better than any others to the detection and measurement of 
minute inequalities in the Earth's motion. 

MDCCCXXXII. R 
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In 1639 . . 


For the Earth. 

. +0",89 . . . 


For Venus. 

. . - 1",28 


1761 . . 


. -0,98 . . . 


. . +1,41 


1769 . . 


. — 1 ,34 . . . 


. . +1,91 


1874 , . 


. . +0,68 . . . 


. . -0,97 


1882 . . 


. . +1,07 • • • 


. . - 1 ,53 



I shall now show the coincidence of the theoretical results with the observa- 
tions that first suggested their necessity. 

From Burckhardt's examination of Maskelyne's observations (Connais- 
sance des Temps, 1816), and from my examination of Mr. Pond's observations 
(Phil. Trans. 1828), it appeared that the mean longitudes of Delambre's tables 
ought to be increased 

in 1783 by 0",25 
in 1801 by 0,08 
in 1821 by 2 ,05 

These observations are all reduced by the same catalogue. The differences of 
the corrections are not proportional to the intervals ; and this is the circum- 
stance that shows the existence of some periodical inequality. 

Now the values of the argument of the long inequality in the epoch are 

for 1783 .... 240° 59' 
for 1801 .... 268 4 



for the middle of 1821 



298 46 



The sines of these angles are — 0,8745, — 0,9994, — 0,8766 ; and hence the 
values of the inequality were 

in 1783 . . . . - 1",80 

in 1801 . . . . -2,06 

in 1821 . . . . -1,81 

If these had been applied in the tables, the corrections given by the observa- 
tions above would have been 



in 1783 . . 


. . 2",05 


in 1801 . . 


. . 2 ,86 


in 1821 . . 


. . 3 ,86 



IN THE MOTIONS OP THE EARTH AND VENUS. 123 

and the differences between these are almost exactly proportional to the times. 
They show that the secular motion ought to be increased by 4",8 (the preces- 
sion being supposed the same as in the application of Maskelyne's catalogue) ; 
and then the application of the inequality investigated in this memoir will give 
correctly the Sun's mean longitude. 

It appears, however, that the inequality in the motion of the perihelion given 
by this investigation, will not account for the anomalies in the place of the 
perihelion given in my paper referred to above. 



Thus terminates one of the most laborious investigations that has yet been 
made in the Planetary Theory. The term in question is a striking instance 
of the importance to which terms, apparently the most insignificant, may some- 
times rise ; and the following remark will show the magnitude of the errors 
which might, under other circumstances, have arisen from the neglect of this 
term. If the perihelia of Venus and the Earth had opposite longitudes, and if 
the line of nodes coincided with the major axes, the excentricities and incli- 
nation having the same values as at present, the coefficient of the inequality 
in the epoch would be 8",9 ; and all the other terms would be important. A 
very small increase of the excentricities and inclination would double or treble 
these inequalities. 

I have avoided any discussion of physical theory, as little can be added 
at present to what has been done by Laplace and others. I may remark, 

ds 
however, that my expression for -r-. differs from that given by Laplace ; and 

that the difference produces no effect in the ultimate result, because Laplace 
nsesfndt where I have usedratf. On this point I have only to state that, by 
adopting the expression which I have used, every formula for the longitude, 
the radius vector, and the velocity in any direction, is exactly the same in 
form for the variable ellipse as for an invariable ellipse (taking the variable 
elements instead of constant ones). If the disturbing force should at any 
instant cease, my value of s for that instant would be the true value of the 
epoch of mean longitude in the orbit which the planet would proceed to 
describe. It is precisely the object of using the method of variation of ele- 

r2 
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merits, to obtain expressions which possess these properties ; and therefore I 
have little doubt that my form will be recognised as more completely in ac- 
cordance with the principles of that method than Laplace's. I should not in 
the present instance have raised a question on this point, but that I conceive 
the method of variation of elements, or some similar method, possessing the 
same advantages of simplicity of application and unlimited accuracy as to the 
order of the disturbing force, will ere long be adopted in the Planetary Theo- 
ries, to the total exclusion of other methods. With this expectation, it appears 
important to adhere closely to the principles of the theory in every formula 
that is derived from it. 

I believe that the paper now presented to the Royal Society contains the 
first* specific improvement in the Solar Tables made in this country since the 
establishment of the Theory of Gravitation. And I have great pleasure in 
reflecting that, after having announced a difficulty detected by observation, I 
have been able to offer an explanation on the grounds of physical theory. 



Postscript. 

In estimating the variation of the elements of the orbit of Venus, the 
change of longitude of perihelion was supposed to be the same as the sidereal 
motion of the perihelion. This is not strictly true ; as the longitude of the 
perihelion, measured as in Art. 4, depends upon the place of the node, and is 
affected therefore by the motion of the node as well as by the motion of the 
perihelion. The amount of the error is however perfectly insignificant. 

G. B. Airy. 
Observatory, Cambridge, 
Nov. 8, 1831. 

* I am not aware that anything has been added to the theory of planetary perturbation, by an 
Englishman, from the publication of Newton's Principia to the communication of Mr. Lubbock's 
Researches. In Maskelyne's tables are two for the perturbations of the Earth produced by Venus 
and Jupiter, calculated Che states) by himself ; but they are utterly useless and erroneous, as they con- 
tain no terms depending on the excentricities. 



